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Azimuth quadrupole component spectra on transverse rapidity y t 
for identified hadrons from Au-Au collisions at ^/snn = 200 GeV 

Thomas A. Trainor 
CENPA 354290, University of Washington, Seattle, WA 98195 
(Dated: November 19, 2008) 

I present the first isolation ol azimuth quadrupole components from published V2(pt) data (called 
elliptic flow) as spectra on transverse rapidity yt for identified pions, kaons and Lambdas/protons 
from minimum-bias Au-Au collisions at 200 GeV. The form of the spectra on y t indicates that 
the three hadron species are emitted from a common boosted source with boost Aj/to ~ 0.6. The 
quadrupole spectra have a Levy form similar to the soft component of the single-particle spectrum, 
but with significantly reduced (~ 0.7x) slope parameters T. Comparison of quadrupole spectra with 
single-particle spectra suggests that the quadrupole component comprises a small fraction (< 5%) 
of the total hadron yield, contradicting the hydrodynamic picture of a thermalized, flowing bulk 
medium. The form of V2(pt) is, within a constant factor, the product of p' t (pt in the boost frame) 
times the ratio of quadrupole spectrum to single-particle spectrum. That ratio in turn implies that 
above 0.5 GeV/c the form of v^ipt) is dominated by the hard component of the single-particle 
spectrum (interpreted as due to minijets). It is therefore unlikely that so-called constituent- quark 
scaling attributed to vi is relevant to soft hadron production mechanisms (e.g., chemical freezeout). 

PACS numbers: 12.38.Qk, 13.87.Fh, 25.75.Aq, 25.75.Bh, 25.75.-q, 25.75.Ld, 25.75.Nq, 25.75.Gz 
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I. INTRODUCTION 

Measurement and interpretation of "elliptic flow" [y%) 
is considered of central importance to the RHIC program 
because it provides the main support for interpretation 
of RHIC collisions as producing sQGP — a thermalized 
strongly- coupled partonic medium with very low viscosity, 
sometimes described as a "perfect liquid" [H, H, Hj] . The 
differential form v 2 (m ) hi particular is a keystone of that 
interpretation d, H, H, 0] . 

In the conventional flow description ua(pt) for identi- 
fied hadrons follows trends predicted by hydrodynamic 
(hydro) models at smaller p t @, H, Q but "saturates" 
at larger p t where parton fragmentation is expected to 
dominate [To[. V2(pt) "scaling" is used to demonstrate 
that "constituent quarks" dominate hadronization. v 2 
and pt divided by constituent quark number n q (i.e., 
2 for mesons, 3 for baryons) appear to be related by 
a universal curve, the inference being that hadrons are 
formed by quark coalescence from a thermalized partonic 
medium [11 1 . Further evidence for collective partonic flow 
is inferred from v 2 data for selected hadrons such as the 
(j> and D mesons and 3 and f2 baryons, where elliptic flow 
generated by hadronic rescattering should be small [l2j . 

The large v 2 at RHIC energies, described in a hydro 
context as elliptic flow, is thus interpreted to imply early 
thermalization of a collective partonic medium resulting 
in large pressure gradients which drive the development 
(for an azimuthally asymmetric system) of the observed 
azimuth eccentricity in hadronic momentum space 0]. 

Logical, technical and interpretational problems have 
emerged for v 2 . Formation of a thermalized partonic 
medium may imply measured v 2 systematics in a hydro 
context, but do v 2 data require the hydro interpretation? 
If more accurate measurements of v 2 or a different asym- 
metry measure are introduced can the hydro interpreta- 



tion be falsified? Is hadron formation from an extended 
QCD field system in some sense a universal characteristic 
of all nuclear collisions at RHIC energies, including N-N 
collisions? Is partonic or hadronic rescattering required 
to produce a system which appears to be thermalized? 

v 2 (pt) as defined is a ratio of two spectra, confusing 
single-particle two-component physics with the physics of 
the azimuth quadrupole (cf. [13| for quadrupole terminol- 
ogy). v 2 data are therefore difficult to interpret directly 
in terms of conventional spectrum analysis. Comparisons 
are typically made indirectly via hydrodynamic models 
whose validity can be questioned, especially because they 
do not model important aspects of single-particle spectra. 

Recent initiatives have shed new light on the azimuth 
quadrupole problem. Single-particle spectra for identi- 
fied hadrons have been accurately separated into soft and 
hard components (longitudinal and transverse fragmen- 
tation) [h|. No evidence for collective radial flow was 
found. An analysis of the algebraic structure of v 2 and al- 
ternative measures reveals that two-particle correlations 
are basic to any v 2 measurement, and r\ dependence of 
2D angular autocorrelations can be used to isolate az- 
imuth quadrupole correlations from "nonflow" [l3|, [l5| . 
Re-examination of the centrality dependence of published 
p t -integrated v 2 data reveals a simple dependence on the 
number of binary collisions, and minijets are identified 
as the dominant source of nonflow [l5| . 

In this analysis spectra associated with the azimuth 
quadrupole are extracted from published v 2 (pt) data and 
plotted on transverse rapidity yt- Those spectra deter- 
mine the quadrupole source boost and relative hadron 
abundances. The main goals of the present analysis are 
to identify the quadrupole component as single-particle 
spectra for each hadron species and to determine the 
abundance of quadrupole hadrons for each species. 

Quadrupole spectra and associated hadron yields could 
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play a critical role in tests of hydro model validity and 
claims of "perfect liquid" in response to open questions. 
For instance, what are the spectrum properties of the 
quadrupole, and what fraction of all produced hadrons 
does the quadrupole component represent? I.e., do al- 
most all particles in a collision participate in elliptic 
flow as widely assumed (a truly collective flowing bulk 
medium), or is the quadrupole component an isolated 
process involving a small fraction of the total system? 
Specific tests of hydro theory are reserved for subsequent 
analysis. 

This paper is arranged as follows. The analysis proce- 
dure is briefly outlined. New methods derived from two- 
particle correlation analysis are reviewed, and conven- 
tional elliptic flow analysis is interpreted in that broader 
context. A sample of 1*2 (pt) data for three hadron species 
is transformed between different plotting formats to il- 
lustrate the requirements for a full analysis and its likely 
outcome. The full analysis includes three steps: 1) spec- 
ify a two-component representation of the single-particle 
spectrum for each hadron species, 2) incorporate the 
kinematics of boosted particle sources, 3) extract and 
interpret azimuth quadrupole spectra. 

Steps 1), 2) and 3) are used to obtain quadrupole 
spectra on transverse rapidity y t from selected v 2 (pt) 
data. Quadrupole spectrum shapes are compared to 
single-particle spectra to search for quadrupole manifes- 
tations therein and to place limits on quadrupole absolute 
yields. Quadrupole and minijet contributions to spectra 
are compared to determine the relations between them, 
and several v 2 (pt) scaling relations and their implications 
for claims of sQGP are examined in the context of this 
analysis. 



II. ANALYSIS DESCRIPTION 

Conventional differential flow measure v 2 (pt) as de- 
fined includes a ratio of two hadron spectra: the spec- 
trum of "flowing particles" (quadrupole component) in 
the numerator and the azimuth-averaged single-particle 
spectrum in the denominator. Although arguments from 
a hydrodynamic context favor ratio vi , it is important to 
examine the quadrupole spectrum (numerator) directly. 

The goal of this analysis is to isolate from existing 
Viipt) data the spectra on y t of hadrons associated with 
the quadrupole component for three hadron species, and 
to compare quadrupole spectra with azimuth-averaged 
single-particle spectra. The analysis should reveal the 
radial boost distribution of the particle source and the 
fractional yields of quadrupole hadrons in a collision. The 
analysis should improve our understanding of the under- 
lying ph ysic al mechanism. E.g., is it hydrodynamic ex- 
pansion jlftl ] or QCD field interactions jl5j? 

A hint of the benefits of this analysis is obtained by 
plotting «2(pt)/pt(lab frame) vs the proper transverse ra- 
pidity y t for each hadron species, as in Sec. HVl To un- 
derstand why that strategy provides qualitatively new in- 



formation the kinematics of boosted thermal sources are 
reviewed in Sec. IVH To isolate the quadrupole spectrum 
(numerator) from the V2(j?t) ratio corresponding single- 
particle spectra (denominator) are presented in Sec. [V] 

In the full analysis Fourier amplitude V 2 {pt) is re- 
covered from V2{pt) by eliminating the single-particle 
spectrum from its denominator. Based on the Cooper- 
Frye description of a thermal source boosted on trans- 
verse rapidity yt a factor p' t (pt in the boost frame) 
is also removed to form an approximate expression for 
quadrupole spectrum p^ivt't ^Vto)- There remains an 
O(l) factor due to an integral approximation and ambi- 
guity between the quadrupole boost Ay t 2 and the abso- 
lute yield of the quadrupole spectrum n c h 2 . Comparisons 
between quadrupole and single-particle spectra, includ- 
ing the hard component (scattered-parton fragments), 
constrain the absolute quadrupole spectrum and yield. 

III. AZIMUTH CORRELATION ANALYSIS 

Two-particle azimuth correlation analysis is outlined, 
and conventional differential (on pt) elliptic flow analysis 
is described in the larger context. Given a pt spectrum 
defined in histogram form with bins of width Spt, the 
symbol Vzipt) indicates the value of v 2 in a pt bin with 
bin multiplicity n Pt . The measured integral quantities for 
each collision event of n particles are vector Fourier co- 
efficients Q m = 53j=i u{m4>i) and scalar power- spectrum 
elements V£ = ES=i"W>) • u(m(pj) [H, The 
same quantities in differential form can be defined as ID 
and 2D histograms respectively on p t bins. 

A. Two-particle correlations on pt 

Two-particle azimuth correlations can be studied with- 
out introducing a reaction or event plane. The basic mea- 
sures of sinusoidal azimuth correlations are the Fourier 
power spectrum elements [l3[ . The 2D p t -integrated 
quadrupole term V 2 2 can be generalized to a pt-differential 
form with unit vectors u(2<f>i). 

n Ptl > n Pt2 

Vi(p tl , P t2) = cos(2[&-&]) (1) 

n Ptl ' n Pt2 

= V 2 (pti) ■ V 2 {pt2), 

where e.g. index p t \ labels a histogram bin of nominal 
width Spt with center at pti containing n pn particles. 
The dot product in the last line defines a mnemonic rep- 
resentation of the i 7^ j double sum. Individual vectors 
Vzipt) ar e not accessible. Diagonal element V 2 (p t ,pt) 
denotes the power spectrum element for a single bin cen- 
tered at pt- V 2 (pti,Pt2) is an additive two-particle cor- 
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relation measure, playing the same role for the azimuth 
quadrupole that total variance ^ Pt:n plays for p t fluctu- 
ations/correlations V 2 2 (pti,Pt 2 ) can describe a two- 
particle distribution on transverse momentum (j>ti,Pt 2 )j 
mass (m i i,m t2 ) or rapidity y t2 )- 



B. Marginal distribution V 2 (pt) vs v 2 (pt) 

Marginal distribution V2{pt) is obtained from the 
asymmetric 2D case that one pt bin is the entire accep- 
tance (including n particles). I.e., V 2 (pt) is obtained by 
integrating V 2 (pti,Pt 2 ) over one Pt axis 



vi( Pt ) 

V 2 (p t ) 
v 2 {2}(p t ) 



7l pj . ,71— 1 

£ 

V 2 {p t )-V 2 
V 2 {pt)-V 2 

v 2 

V 2 {p t )/n Pt 
V 2 (pt)-V 2 



'■pt 



Vn 



The last line defines elliptic flow measure v 2 {2}(p t ) in 
terms of two-particle correlations. In general, V 2 (pt) ^ 
\ZV 2 (PtiPt)- V 2 (pt) is an element of the marginal distri- 
bution, whereas V 2 2 (pt,Pt) refers to a single diagonal bin 
on (pti,Pt 2 )- The two are related by a covariance. 



C. Conventional event-plane method 

Conventional v 2 analysis is motivated in the context of 
an event or reaction plane, but analysis results do not de- 
pend on a reaction plane per se. v 2 measures the m = 2 
Fourier component of any two-particle azimuth correla- 
tions present in collision products, including jet correla- 
tions. An "event plane" can arise from any event- wise 
azimuth structure (including minijets), and the "event- 
plane resolution" may not relate to a true reaction plane. 

Pi-diffcrcntial elliptic flow analysis at mid-rapidity is 
based on a ID Fourier decomposition on azimuth of the 
?7-averaged 3D density. The Fourier series is defined in 
terms of reaction-plane angle ^ r 



p(Pt, 4>) = \ 1 + 2 ^2 y ™{pt) cos(m[0 - 



2tt 




(3) 



where Vo(pt) /2tt = po(pt) is the 3D single-particle 
pt spectrum (averaged over 2tt azimuth and one unit 
of pseudorapidity about n = 0) described by a two- 
component spectrum model jl4j [. Fourier amplitude ra- 
tios v m { Pt ) = V m {p t )/V (p t ) = (cos{n#-* r ]}(p t )) EH- 
V m could represent multiple physical contributions, in- 
cluding minijets as well as various "flow" sources. Eq. ([3]) 



is not a conventional Fourier series because common ele- 
ment po divides each term, thereby coupling all v m . The 
equation is nonphysical, since \Ev is not known a priori, 
and the V m are therefore not measurable by inversion. 

Within the flow model description \t r must be esti- 
mated from a subset of the collision products. The esti- 
mate is called the event-plane angle V&m, and Eq ([3]) is 
rewritten in terms of unit vectors uirncf)) as 



p{pt,4>) = £ Qm(Pt)-u(r 



(4) 



m—~oo 



= ^ i 1 + 2 E 9m(jH) cos (mfe - * m ]) \ , 



m— 1 



(2) with true Fourier coefficients Q m {pt) = Sjep t ^( m 0j) = 
Qm(pt) u(mty m [p t ]) and Fourier amplitude ratios 
Imipt) = Qm{pt)/Qo(Pt)- The Q m are conventionally 
interpreted by assuming that azimuth structure is 
hydrodynamic in origin (various flows) relating to 
the reaction plane. However, the Q m may contain 
substantial "nonflow" contributions dominated by the 
Fourier coefficients of the same-side minijet peak (jet 
cone) 15]. The inferred "event-plane angle" ^ m (actu- 
ally the Fourier phase angle) may be poorly correlated 
or uncorrelated with the actual A-A reaction plane. 

The differential amplitude ratio q 2 (pt) can be obtained 
by inverting the Fourier series 



Q 2 (Pt 



{u(2ct> iept ) ■ u{2y 2 \p t })) 

Q 2 (pt) Q 2 (pt) _ Qiipt) 



(5) 



i(pt) 



with n Pt — Qo(pt)- However, according to standard flow- 
analysis methods Q 2 {pt)/Q 2 {pt) — > Q 2 /Q 2 = u(2^ 2 ) 
which determines the (global?) m = 2 event-plane angle, 
and "autocorrelations" (self pairs) must be eliminated 
from the dot product [191 ] . For each particle i in a pt bin 
a complementary vector Q 2 — > Qn is formed by omit- 
ting the i th particle from the Q 2 sum over j, q 2 (pt) then 
becomes conventional elliptic flow measure v 2o b s {pt ) 

v 2obs (pt) = (u(2(j) iept ) -u(2*2i)) (6) 

= I QM Q 2i \ 
\ ript Q 2l I 

M V 2 (pt)-V 2 
n pt (Q 2l ) 

where the V 2 dot product defined in Eq. |2j) represents 
the double sum with j =/= i. v 2o bs is then divided by the 
"event-plane resolution" (cos(2[4 , 2 — ^r])) to obtain 



v 2 {EP}(p t ) 



(u(2<f> iept )-u{2t> 2l )) 
(cos(2[* 2 - * P ])) 

V 2 /(Q 2i 



(7) 



M2}(Pt) 



(cos(2[* 2 -* r ])) : 
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which gives the exact relation between V2{EP} and i>oj 2| 
for the first time in terms of the 0(1) second factor [L3| . 
The difference between {EP} and {2} results from a mis- 
conception about the V2obs numerator leading to intro- 
duction of (Q2i) ~ Q 2 in the denominator of Eq. ([6]) in 
place of V 2 as in Eq. ([2]) (last line) [13, [ll]. The extrane- 
ous "event-plane resolution" ~ V 2 / Q2 is then introduced 
to correct v 2obs (p t ). 

The event- plane method is also described in terms of 
"subevents," [2(|. A correlation quantity is defined by 

^ n Pt n — l 



"Pt 



i&Pt 



V 2 (pt)-V 2 



"■Pt 

and normalization is obtained from 

Q 2a -Q2b = cos(2[4-0,]), (9) 

where a, b denote two equivalent and disjoint partition el- 
ements (subevents) covering a detector acceptance (n a — 
rib ~ n/2). If the disjoint partition elements are perfectly 
correlated, and there are no nonflow contributions, then 
Q 2a ■ Q2b = V 2a ■ V 2b ~ V 2 2 /A [TEEl and 



(u(2<t>, 



£pt , 



Q2 



2yQ 2a -Q 



21, 



YM ■ V2 

n Pt V 2 



v 2 {2}( Pt ), (10) 



explaining Eq. (5) of [20|. If Eq. (JTOJ) is multiplied top 
and bottom by 1/Q 2 the EP method Eq. is approx- 
imated. The first relation in Eq. (|10p is approximate 
because of uncertain physical implications of the defini- 
tion of (cut system for) partition elements a and b. The 
definition of the (a, b) partition may reduce nonflow con- 
tributions to Q 2a ■ Q 2 b compared to V 2 (pt) ■ V 2 , leading 
to an undetermined systematic error in ratio v 2 (pt). 

This section demonstrates that v 2 {EP}(p t ) 
v 2 {2}(pt) approximates a generic two-particle correlation 
analysis on azimuth, although motivating language and 
symbols (e.g., "flow vector" Q 2 ) imply that the event- 
plane method necessarily relates to "collective" (hydro- 
dynamic) phenomena. Multiplicities elements of 
the histogrammed single-particle pt spectrum. The pt 
spectrum in the v 2 (pt) denominator obscures interpreta- 
tion and comparisons to theory, as shown in this analysis. 



IV. v 2 (p t ) MEASUREMENTS 



pattern of flows. Larger-p t hadrons from parton frag- 
mentation (possibly by coalescence of "shower" partons) 
reveal modification of fragmentation by the medium. 
Intermediate-pt hadrons may result from recombination 
of "thermal" and "shower" partons [ill l2ll I22I I23I]. 

In the present analysis qualitative conceptual issues are 
of central importance. A simple and accurate data sam- 
ple including both p t and mass dependence is used to 
demonstrate the algebraic structure of v 2 (pt) and the 
basic properties of quadrupole spectra. Notable theory 
examples are included to explore the general relation of 
hydro theory to the quadrupole component in different 
manifestations. v 2 {pt) data for pions, kaons and Lamb- 
das are related to single-particle spectra for pions, kaons 
(interpolated) and protons. Proton and Lambda spec- 
trum shapes are assumed equivalent for this analysis. 

Fig. [1] (left panel) shows data from v 2 {pt) analysis of 
identified mesons (pions, kaons) and baryons (Lambdas) 
from minimum-bias 200 GeV Au-Au collisions [2J, |25j . 
The mass trend at smaller p t (massive hadrons have 
smaller v 2 ) is commonly interpreted to imply collective 
flow (hydrodynamics). At larger p t v 2 data are said to 
"saturate," following a nearly constant trend beyond 4 
GeV/c [nj. 

Hydrodynamic models provide a semi-quantitative de- 
scription at smaller p t but fail at larger p t (hydro mod- 
els overpredict v 2 at larger p t ) [2(|. The dotted curves 
in each panel are viscous hydro predictions with zero- 
viscosity limit for pion v 2 (jpt) (A from 1], B from [27l]). 
The systematics of viscous hydro predictions compared 
to data are interpreted as evidence for a fluid medium 
with very small viscosity ( "perfect liquid" ) P, [1, Q . 




p t (lab) (GeV/c) 



y t (7t,K,p) 



FIG. 1: Left panel: V2( pt) data for three hadron species plot- 
ted in the usual format [24, 25]. Right panel: The same V2(pt) 
data divided by pt in the lab frame suggest a universality on 
proper transverse rapidity for each hadron species — in partic- 
ular the correspondence of data near yt = 1. Dotted curves 
A and B in each panel are viscous hydro predictions from 1] 
and [13] respectively. The three curves through data are de- 
rived from this analysis. 



v 2 data are described in the context of a thermal- 
ized, collectively-flowing bulk partonic medium probed 
by flow measurements and energetic scattered partons. 
Smaller-j^ hadrons emerging from the bulk medium (pos- 
sibly by coalescence of "thermal" partons) exhibit a 



Fig. [T] (right panel) shows the same data in the form 
v 2 {pt) / Pt{\sh frame) plotted vs yt{it, K,p) (proper rapid- 
ity for each hadron species), where transverse rapidity 
yt = log{(m t +pt)/mo}. The first Lambda point at 
yt ~ 0.57 (not visible) is slightly negative but consistent 
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with zero. The simple transformation, revealing peaked 
distributions with similar amplitudes and common left 
edges for the three hadron species, suggests that more 
information can be extracted from existing vi data with 
a generalized analysis method. 

Data distributions in the right panel taken together 
imply that the three hadron species are emitted from a 
common moving (boosted) source, as demonstrated be- 
low. The three curves from this analysis passing through 
data in each panel are based on that hypothesis. The 
relevant model parameters are summarized in the panel. 
The dotted hydro curves in the right panel p], [27j devi- 
ate significantly from the pion data trend. However, the 
same curves compared to pion data in the left panel have 
been cited to imply small medium viscosity and forma- 
tion of a "perfect liquid" at RHIC. The relation of the 
hydro curves to data is discussed in Sec. IXII 



V. SINGLE-PARTICLE SPECTRA 

The first step of this analysis is to obtain the single- 
particle azimuth-averaged 3D spectrum po = Qo(yt)/2n 
(denominator of i>2[pt]) for three hadron species from Au- 
Au collisions at 200 GeV which provides a context for 
the quadrupole component (numerator of i>2[pt])- Single- 
particle spectra may include a quadrupole contribution 
which can be used to estimate the absolute quadrupole 
yield. Two-component spectrum analyses of p-p and Au- 
Au spectra reported in [3, S3 are use d to construct 
minimum-bias spectra compatible with the data in Fig. [1] 



A. Spectrum notation 

The 3D single-hadron density on momentum averaged 
over one unit of pseudorapidity rj about mid-rapidity is 

1 d 3 n 

Po(xt,(p) = — (11) 
Xt dxt drj d<p 

^ ^ _ 1 d 2 n 

2irxt dxt drj ' 

where the second line is averaged over azimuth. Trans- 
verse measure x± is pt, rat or y t . Transformations be- 
tween densities require Jacobians dyt/dp t — l/ra t and 
dy t /dm t = 1/pt- For reference, p^N = d i^nn /dr/dtfi ~ 
2.5/27T is the rj- and ^-averaged, pj-integrated 2D hadron 
density at mid-rapidity for 200 GeV NSD N-N collisions. 

It is sometimes useful to plot all hadron species on 
pion rapidity denoted by yt-n ~ ln(2p t /m Tr ) or (for plot 
axes) 2/t(7r). yt-n is then simply a logarithmic measure of 
Pt providing better visual access to spectrum structure. 
When relativistic transformations (boosts) are important 
the proper y t for each hadron species should be used, 
denoted by variable yt with no qualification and plot axis 
labels y t (ir,K,p). 



Comparison of results on transverse variables pt, rat 
and yt, as in this study, is essential to distinguish dif- 
ferent physical mechanisms. For thermal spectra rn t is 
preferred. For boosted systems proper y t for each hadron 
species is preferred. For parton fragmentation (minijets) 
Pt would reflect the common underlying parton spectrum, 
but ytTT (~ ln(2p t /m,r)) provides better visual access to 
structure. Analysis of spectra on a single plotting vari- 
able may confuse several dynamical mechanisms. 



B. Glauber model and multiplicity definitions 

The Glauber model of A-A collisions defines several 
A-A geometry parameters (2^. For A-A impact param- 
eter b, np ar t/2 is the corresponding average number of 
participant nucleon pairs and nu n is the average num- 
ber of N-N binary collisions (for a specified scattering 
process). Some hadron production processes are propor- 
tional to n part /2 (soft), and some are proportional to nun 
(hard). The combination comprises the two-component 
model of hadron production, which describes N-N col- 
lisions well [28| and serves as a reference in A-A colli- 
sions [14| . v = 2nbinl n parti the mean participant path- 
length in number of encountered nucleons, is a geometry 
parameter used to measure A-A centrality. 

n c h is the total charged-particle multiplicity in one unit 
of rj at mid-rapidity. The total multiplicity associated 
with the quadrupole component is n c h2- The quadrupole 
multiplicity associated with hadron species X is n c hX2- 
Quadrupole multiplicities should not be confused with 
quadrupole Levy distribution shape parameter ri2 or nx2- 

Ambiguities in the normalizations of measured spectra 
are noted in [141] . specifically the centrality dependence 
of integrated n c h compared between experiments. The 
present analysis concerns minimum-bias V2(pt) data for 
which an average over centrality is implicit. The asso- 
ciated normalization uncertainty in the averaged single- 
particle spectra is about 20%. However, normalization 
uncertainty is not relevant to the present analysis which 
refers only to relative spectrum shapes. 



C. Two-component spectrum model 

The two-component (soft+hard) model of hadron spec- 
tra provides a compact and accurate description of p-p 
and Au- Au collisions [3, HH . The soft component is in- 
terpreted as longitudinal participant-nucleon fragmenta- 
tion. The hard component at mid-rapidity is interpreted 
as minimum-bias large-angle scattered parton fragmen- 
tation (minijets), which can also be described as hadrons 
emitted from a radially-boosted source. The open ques- 
tion for any observed boost phenomenon is what physical 
mechanism produced the boost. 

The two-component models of pion, kaon and proton 
spectra (3D densities per participant pair) at 200 GeV 
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are summarized by 



-POrr = f {Son +0.012 VTAA-K #0tt} (12) 



Tlpart 

2 

Tlpart 

2 



-POK 



POp 



1.012 
0-09 PNN 

1.16 
0.06 pnn 

1.12 



{S^if + 0.16 VT A AK H 0K } 
{S 0p + 0.l2isr AAp H Qp }, 



with the differential form of pox(yt) defined in Eq. (fTTj) . 
Unit-integral model functions Sox(yt) and H x(yt) and 
hard-component ratios taax^Vu v) for pions and protons 
are defined in 
from the N-N - 



14] 



The r aax represent all deviations 
Glauber two-component linear reference. 
Kaon model functions were estimated by interpolation for 
this analysis. 




FIG. 2: (Color online) Summary of pion (left) and proton 
(right) per-participant-pair single-particle spectra from Au- 
Au collisions at 200 GeV and five centralities [l4|. Hnn is the 
hard component (minimum-bias transverse parton fragmen- 
tation) and Snn is the soft component (longitudinal nucleon 
fragmentation), both inferred for N-N collisions. The solid 
points in the left panel represent the NSD p-p spectrum [28j |. 



Model spectra describing pion and proton data are 
summarized in Fig. [2] Reference soft components So 
(unit-normal distributions not shown) are Levy distri- 
butions on m t transformed to pion yt- The transforma- 
tion strategy is discussed in [1J] . Reference hard compo- 
nents Hq (also unit-normal distributions not shown) are 



Gaussians on y t with exponential tails cx exp(- 



Vt) 



representing expected QCD power law p~[ nh required by 
data above p t ~ 6 GeV/c (y tir ~ 4.5). Distributions 
Snn and Hnn have the same forms but integrate to 
hard- and soft-component hadron numbers n s and Uh 
with n s + nh = rich for N-N collisions [28|. By hypoth- 
esis, the soft component for Au-Au collisions remains 
fixed at the N-N reference. Deviations of hard compo- 
nent ti aa from its N-N reference are measured by ratio 
taa = Haa/ Hnn- The model functions describe the 
shapes of the data spectra at the few-percent level over 
the yt interval relevant to this analysis. The quality of 
the description is indicated by the relative residuals in 
Fig.H 
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FIG. 3: Relative residuals (data — model) / model for pions 
(left) and protons (right) from five centralities of Au-Au col- 
lisions at y / ijvjv = 200 GeV. The two-component spectrum 
model with modification factors taa inferred in [T3] describes 
data to the statistical limits from 0.5 to 10 GeV/c (y t £ [2, 5]). 



D. The quadrupole spectrum component 

Interpretation of the azimuth quadrupole spectrum 
raises a significant question: Does the azimuth 
quadrupole (v% data) represent modulation of a spectrum 
component existing in N-N collisions (e.g., the soft com- 
ponent), or does a new radially-boosted net source of 
hadrons modulated on azimuth emerge in A-A collisions? 
Does that component extrapolate back to N-N collisions? 

The ?7-averaged three-component 3D spectrum on 
(xt , </>) for xt = m t or y t can be expressed as 



p{mt,4>) 

p(yt,4>) 



Po(m t ;T ) 
Po(yf,Po) 



f p2(mt;T 2 ,Pt[4>]) 
- p 2 (yt;p2, Ay t [0]), 



(13) 



where pi is a possible quadrupole (third) component from 
a radially-boosted source. Parameters (3t(4>) and Ay t (4>) 
represent a conjectured azimuth-dependent radial boost 
of the third component. The first term po(yt!A i o) is the 
two-component spectrum from [3| • Quadrupole term pi 
may represent a new particle source, a modification of the 
N-N soft component, of the hard component, or an in- 
teraction between them. To clarify we must estimate the 
shape and absolute magnitude of the quadrupole spec- 
trum component from v^ijpt) data and compare them 
with measured azimuth-averaged yt spectra. 



VI. BOOSTED HADRON SOURCES 

The second step of this analysis is to define the kine- 
matics of nearly-thermal hadron spectra from mo ving 
(boosted) sources, essentially the blast-wave model |3(X 
l3l| related to the Cooper- Frye description of moving (ex- 
panding) particle sources [32j. I consider only monopole 
and azimuth-quadrupole pt and yt spectrum components. 
For simplicity "thermal" spectra are described in the 
boosted frame by Maxwell-Boltzmann exponentials on 
m t . The description can be generalized to Levy distribu- 
tions on mt for accurate descriptions of data. The intent 
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is to provide a general description of hadron production 
from a source including (but not restricted to) a radially- 
boosted component with azimuth variation. 

A. Radial boost kinematics 

The four-momentum components of a boosted source 
are first related to transverse rapidity yt- The boost 
distribution is assumed to be a single value for sim- 
plicity. The particle four-momentum components are 
n%t = mocosh(y t ) and pt — TOosinh(y t ). The source 
four- velocity (boost) components are j t — cosh(Ay t ) and 
7t/?t = sinh(Ay t ), with /3 t = tanh(Ayt). Boost-frame 
variables are defined in terms of lab-frame variables by 

m' t = to cosh(y t - Ay t ) = j t (mt - Pt Pt) (14) 

= m t 7 t {l - tanh(y t ) tanh(Ay t )} 

p' t ee to sinh(y t - Ay t ) = y t (p t - /3 t m t ) 

= m t 7 t {tanh(y t ) - tanh(Ay 4 )} 

with p' t denoted p t (boost) in figures. 

Fig. 3] (left panel) relates p' t — > pt (boost) to p t — > 
Pt(lab). The main source of the mass trend of V2(pt) 
at small pt, interpreted as "hydro" behavior, is a simple 
kinematic effect as seen at lower left. The mass system- 
atics hold for any boosted nearly-thermal hadron source 
independent of boost mechanism (i.e., hydrodynamics is 
not required). The intercepts (p' t = 0) of the three curves, 
given by pto = too sinh(Ay t ), are important for discus- 
sion of the hydro interpretation of V2{pt)- 




p t (lab) (GeV/c) y t (n,K,p) 



FIG. 4: Left panel: p' t (pt in the boost frame) vs p t in the lab 
frame. Factor 7t(l — fit) in the denominator insures that the 
combination — > pt for large pt. Right panel: The same ratio is 
plotted on proper yt for each hadron species, demonstrating 
a fundamental relationship applicable to any hadron species. 

Fig. |4] (right panel) relates p' t to transverse rapidity 
yt (?r, K, p) and illustrates one reason why plots on y t are 
a major improvement over p t or toj. Normalized p' t /pt 

Pt = 1 - ft/ tanh(g/ t ) (15] 

Ptlt(l-Pt) " 1 - Pt 1 } 

increases from zero at monopole boost Ay t o and follows a 
universal curve on y t to unit value for any hadron species. 



Thus, normalized p' t goes asymptotically to pt for large 
Pt (or y t ) independent of boost. The form in Fig. 2] (right 
panel) is important for interpreting v^ipt) data as spectra. 
In the present study we find Ay f o ~ 0.6 ~ 7t (1 — fit), 
common to three hadron species. 



B. Radially-boosted thermal spectra 

The simple blast-wave model invoked here assumes 
longitudinal-boost-invariant normal emission from an ex- 
panding cylinder, eliminating the need for Bessel func- 
tions K\ and Iq arising from integrals over y z and (f> [l6| . 
Slope parameter T for mt spectra and thermal parameter 
(j, = m /T for y t spectra are defined. Boosted spectra on 
yt and mt (32} are 

p{y t \n, Ay t ) = A yt exp{-^i [cosh(y 4 - Ay t ) - 1]} (16) 
p(mt]T,Pt) = A mt exp{-[7t (m t - fi t p t ) - m]/T}, 

providing a simplified description of "thermal" radiation 
from a radially-boosted cylindrical source. Applications 
require a specific radial boost model Ay t (r, <j>). 



C. Radial boost models 

In a nuclear collision there are (at least) two pos- 
sibilities for the radial boost model: 1) a monolithic, 
thermalized, collectively-flowing hadron source ("bulk 
medium" ) with complex transverse flow distribution 
dominated by monopole (radial flow or Hubble expan- 
sion) and quadrupole (elliptic flow) azimuth compo- 
nents [3l| : and 2) multiple hadron sources, some with 
azimuth- modulated transverse boost. Hadrons may 
emerge from a radially-fixed source (soft component), 
from parton fragmentation (hard component), and pos- 
sibly from a source with radial boost varying smoothly 
on azimuth, including monopole and quadrupole com- 
ponents. Case 2 is assumed for this analysis, but both 
possibilities are reconsidered in light of analysis results. 

A radial boost with monopole and quadrupole compo- 
nents is described by 

Ay t (<t>) = Ay t0 + Ay t2 cos(2A0 r ) (17) 
ft fa) = tanh(Ay t [$) 

- fto + Aa cos(2AeV), 

with Ay t 2 < Ay t o for positive- definite boost. The conven- 
tion Acf> r ee (j)— is adopted for more compact notation. 
Monopole boost component Ay t o is easy to extract from 
data, but quadrupole component Ay t 2 is less accessible. 
Ay 4 o could be interpreted as a "radial flow" but may 
apply to only a small fraction of produced hadrons. The 
quadrupole boost magnitude should reflect the eccentric- 
ity e of the A-A collision geometry. 
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VII. AZIMUTH QUADRUPOLE COMPONENT 

The third step of this analysis is to relate the azimuth 
quadrupole spectrum component algebraically to experi- 
mental i>2 data. I assume that 1) the quadrupole compo- 
nent arises from a hadron source with azimuth-dependent 
radial boost distribution Ay t ((f>), 2) the quadrupole 
source may produce only a small fraction of the hadrons 
in a collision, and 3) the quadrupole spectrum may ap- 
pear to be thermal in the boost frame and may be inde- 
pendent of the soft and hard spectrum components. 

A. Quadrupole-component model 

Given those assumptions the ^-averaged 3D spec- 
trum at midrapidity for hadrons associated with the 
quadrupole component is modeled by 

P2(yt,4>)=A 2yt cxp{-^ 2 [cosh(yt - Ay f (») - 1]} (18) 
p 2 (mt, <f>)=A 2mt cxp{- {-f t {(j>)[m t - /3t(4>)Pt] ~ m o) /T 2 }, 

where a Maxwell-Boltzmann (M-B) distribution for a 
locally-thermal source is assumed for simplicity, and 
jj, 2 = TO0/T2 for the quadrupole. The procedure below 
may be applied to a more general function such as a Levy 
distribution. Ay t (4>) defined in Eq. (fT7|) represents fixed 
monopole and quadrupole boost components. 




FIG. 5: Left panel: (Color online) The quadrupole component 
modeled by a thermal spectrum boosted by Ayt (4>) containing 
monopole and quadrupole terms: monopole boost Ay t o = 1 
and quadrupole boost amplitude Ayt2 = 0.5. Right panel: 
projection of the left panel onto yt revealing the edge struc- 
ture. 

Fig. [5] (left panel) illustrates the form of p 2 (yt,<?!>) 
relative to reference or reaction-plane angle , 5 r , with 
Ayto = 1 an d Ay t2 = 0.5. Fig. [5] (right panel) shows 
the projection of P2(Vu4') onto yt- The half-maximum 
point of the left edge of the quadrupole spectrum is at 
monopole boost Ayto- The projection suggests that the 
shape of the left edge might reveal quadrupole boost am- 
plitude Ay t 2 if resolved accurately. However, the small- 
yt region is experimentally difficult. The left edge is also 
affected by variations of Ay t o within a centrality bin. Ac- 
curate data would be needed to determine edge details. 



B. Quadrupole Fourier amplitude 

Fig. [6] (left panels) shows unit-amplitude p 2 (yt,</>) f° r 
Ayto = 0.5 and quadrupole boost amplitude Ay t2 = 
0.125, 0.250. We can now obtain the relation be- 
tween inferred quadrupole Fourier amplitude V 2 (yt) and 
quadrupole spectrum p 2 (y t ; Ayto). The Fourier ampli- 
tude is defined by 

l/ 2 (yt;Ayto,Ay t2 ) = f d^p 2 (y t , <t>) cos(2A0 r ) (19) 

J — 7T 

assuming that reaction-plane angle \& r is known and 
P2(yt,4') represents all m = 2 azimuth dependence in 
the single-particle spectrum. 




FIG. 6: Left panels: (Color online) The quadrupole compo- 
nent modeled as in Fig. [5] for Ayto = 0.5 and two values 
of quadrupole boost amplitude Ayt2- Right panels: Corre- 
sponding Fourier amplitudes V2 (yt) from Eq. (|19[1 normalized 
by quadrupole boost amplitudes Ayt2. Widths of the nega- 
tive regions are 2Ayt2- The dash-dot curve is the quadrupole 
spectrum p 2 (yt; Ay w )- 

Fig. [6] (right panels) shows the integral in Eq. (fl9|) 
times T 2 /27rAy t2 as a histogram (points and thin solid 
curve). The Fourier amplitudes peak at yt ~ 1.5 and 
fall toward zero at yt = Ay t o. Similar amplitudes at 
the peak confirm that the integral is cx Ay t2 . A negative 
undershoot centered at Ayto (vertical dash-dot line) with 
width ~ 2Ay t2 appears because the phase of the sinusoid 
changes by tt in traversing from one side of the mode of p 2 
on y t (y t ~ Ay t0 + Ay t2 ) to the other. Negative values 
of d 2 do not require collective flow of a medium, only 
boost of the quadrupole component. The other curves 
are described below. 
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C. Factoring V 2 (yt) 



Combining factors we obtain 



Eq. ([19]) can be factored to isolate the underlying 
quadrupole spectrum p 2 (yt ! Ay*o)j the subject of this pa- 
per. Invoking the Ay t (<fi) model defined above and re- 
ferring back to Eq. (|T4)) I expand the cosh term in the 
boosted M-B distribution of Eq. (fig)) as 



cosh(y t - Ay* [</>]) -1 = cosh(y t - Ay t0 ) - 1 + (20) 
cosh(y t - Ay t0 ){cosh(Ay t2 cos[2A</V]) - 1} + 
sinh(y t - Aj/ t o) sinh(Aj/ t2 cos[2A</> r ]) 
The three terms correspond to three factors of p 2 (yt, 4 1 )- 
p 2 {yt,4>) = M exp{-(m;-m )/T 2 } x (21) 
exp{m' t [cosh(Ay t2 cos[2A0 r ]) — 1]/T 2 } x 
exp{p' t sinh(Ay t2 cos[2A0 r ])/T 2 } 
= p 2 (yt;Ay t o) x Fi(y t ,<f>; At/ t o, Ay t2 ) x 
F 2 (yt,<j>; Ay to, Ay t2 ). 

The last line defines azimuth-dependent factors F\ (yt , <p) 
and F 2 (yt , 4>) in terms of monopole and quadrupole 
components of the radial boost. The objective is 
quadrupole spectrum component p 2 (y t ; Ay t o) emitted 
from the boosted particle source as one factor of mea- 
sured Fourier amplitude V 2 (p t ) inferred from v 2 (p t ) data. 

The quadrupole boost dependence is contained in inte- 
S ral 2¥ I-„ dtpFxiyt, 0) F 2 (y u cf>) cos(2A0 r ). The leading 
azimuth dependence of F\ — 1 is cos 2 (2A0 r ) which does 
not contribute appreciably to the integral, so F\ ~ 1 
is a good approximation. Fig. [7] (left panel) shows 
T 2 (F 2 - l)/p' t cos(2A(j) r ) ~ Ay t2 = 0.2. Additional az- 
imuth structure due to higher-order terms in the expo- 
nential is substantial at larger y t . The full integral over 




2 3 

y t (ji,K,p) 



FIG. 7: Left panel: (Color online) The structure of factor 
F2 in Eq. (|2ip . dominated by Ayt2 at smaller y t . Right 
panel: The structure of O(l) factor f(yt\ Ay t o, Ay t 2) defined 
in Eq. Qg). 

factors Fx and F 2 is 

i- JjtFWFM) cos(2A<f> r )^ Ayt2 ^ iyt \ (22) 

defining f(y t ; Ay t0 , Ay t2 ) as an 0(1) correction factor 
plotted in Fig. [7] (right panel) for a particular combi- 
nation (Ay tQ , Ay t2 ). f(y t ) is closer to 1 the smaller is 
Ay t2 /Ay m . 



V 2 (yt;Ayto,Ay t2 ) 
2tt 



Ay t2 p 2 (y t ;A 2/i0 ).(23) 



In Fig. [6] (right panels) the dashed curves through points 
represent p' t /2 ■ f(y t ) p 2 (y t ; Ay t0 ), which agree well with 
the direct integrals (points) except in the region of 
negative values near Ay t o where accurate comparisons 
with data are not possible. The dash-dot curves rep- 
resent p't/2 ■ p 2 {yu Ayto) which does not include factor 
f(y t ; Ay t2 , Ayto). Small deviations from the exact inte- 
gral (dashed curves and points) due to omission of factor 
f(yt) depend on ratio Ay t2 /Ay t o as noted. 

Because V 2 /27r = pov 2 , v 2 {pt) can be expressed as 



v 2 (pt) 



P'i 



f (Vt ) Ay t 2 p 2 (y t ; Ay t0 , T 2 , n 2 ) 



2T 2 



po(yf,T ,n ) 



(24) 



Given Eq. (j2"4")l we can reconstruct at least the shape of 
quadrupole spectrum p 2 (y t ; Ay t o) from measured v 2 (p t ) 
data. Although the derivation is based on an exponential 
form for p 2 , the procedure can be applied to the more 
general form of a Levy distribution within the limited pt 
range relevant to quadrupole and soft components (< 2 
GeV/c in the boosted frame). 



D. Obtaining P2{yt\ Ay t o) from measured V2 data 

The quadrupole spectrum is best related to measured 
quantities with the equation 



po(yt) 



V2(yt) 



Pt 



Pi 



Ptlt(l-Pt)} I 2T 2 
fiVu Ayto, Ay t2 ) Ay t2 p 2 (y t ; Ay t0 ) 



x(25) 



Quantities on the LHS are measured experimentally. 
p 2 (yt', Ayto) on the RHS is the sought-after quadrupole 
spectrum. The common monopole boost Ay t o and T 2 for 
each hadron species can be estimated accurately from the 
p 2 (yt', Ayto) spectrum common left edge and shape. As 
shown in Fig. 0] (right panel) p' t /{ptjt{i — Pt)} is deter- 
mined only by Ay t o and deviates from unity only near 
that point. The numerator of the second factor is also 
determined by Ay to- Thus, all factors on the first line of 
the RHS and the shape of p 2 (yt] Ayto) are determined by 
data on the LHS. 

In the second line of Eq. (|25p RHS there is an ambigu- 
ity in the product of Ay t2 and the p 2 (yt, Ayto) amplitude. 
Comparison of the inferred quadrupole p 2 (yu Ay to) spec- 
trum shape, especially the leading edge of the spectrum, 
with measured azimuth-averaged spectrum po for each 
hadron species may place a lower limit on Ay t2 . The 
upper limit Ay t2 < Ay t o assumes positive-definite trans- 
verse boosts. The two limits establish an allowed range 
for quadrupole spectrum integral n c ^ 2 . Ay t2 should be 
common to all hadron species emitted from a boosted 
hadron source, possibly reducing systematic uncertainty. 
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E. Analysis summary 

To summarize, given measurements of v^ipt) the 
quadrupole spectrum is determined with minimal sys- 
tematic uncertainty by the following steps for each 
hadron species 

1. Parameterize single-particle spectrum po(ytn)', ob- 
tain the value of po(ytTi) for each V2 datum 

2. Calculate and plot Poivt-w) V2(j>t)/pt(}ab frame) 

3. Model P2 (y* , Aj/to) by a boosted Levy distribution 

4. Use the model to plot the product of the first and 
last factors of Eq. ([25)) RHS on y t7r 

5. Compare 4) with 2) to determine monopole boost 
Ay t o and temperature Ti plus Levy n-i 

6. Obtain product f(y t ) Ay t 2 P2(yt, Ayto) from 
Eq. (HU) 

7. Compare the inferred pvivti Ayto) shape with 
single-particle spectra to obtain an upper bound 
on the p2(yt,Ayto) amplitude <J4> lower bound on 
Ay t 2 

8. Iterate Ay t2 to optimize f(y t ; Ay t0 , Ay t2 ) 

9. Obtain corrected pi{yt, Ayto) 

Step 7) should include comparison of the approximate 
centrality variation of p2{yt,Ay t o) with the measured 
centrality variations of the single-particle spectrum com- 
ponents on y t TT [14] to tighten constraints. 

VIII. QUADRUPOLE OBTAINED FROM DATA 

The analysis procedure can be illustrated with the data 
shown in Fig. [1] following the steps described in the pre- 
vious section. The procedure is intended to minimize 
model assumptions and systematic errors. Step 1, defin- 
ing single-particle spectrum parameterizations, was de- 
scribed in Sec. |Vl 

A. Forming the LHS of Eq. ([25]) 

Steps 1 and 2 of the analysis produce the LHS of 
Eq. (j2"5")) from V2(pt) data. Fig. [1] (left panel) shows 
the original V2 data in a form which provides little di- 
rect indication of the underlying physics. The left-most 
measured Lambda point (not plotted) is negative, the 
reason now apparent from the discussion in Sec. IVIIBI 
Fig. [T] (right panel) shows V2(pt)/pt which hints at a sim- 
ple boost phenomenon. The common left edge provides 
an initial estimate of Ay t Q. 

Fig. [5] illustrates how to match po(yt7r) parameteriza- 
tions to V2(pt)/pt data (step 1). Single-particle spectra 
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FIG. 8: Formation of quadrupole spectra from V2 (pt) data and 
single-particle spectra - measured quantities combined with 
full-spectrum two-component parameterizations. The open 
symbols are the values of po(yt) used for the conversion. The 
solid symbols are the resulting approximations to quadrupole 
spectra. The dashed curves are from the present analysis. 
The solid curves result from removing a p' t /pt factor. The 
dash-dotted curves are hard-component models from [T3 |. 



in the form 2/n par t ■ Po(ytir) for three hadron species 
(protons and pions from [14j . kaons interpolated) are 
given by the dotted curves [defined in Eq. (fT2"|) and dis- 
cussed in Sec. IV Cj . The open symbols show the spe- 
cific values of po for each V2 datum and hadron species. 
The solid symbols show the corresponding values of 
2/n par t ■ Po{ytir) v 2 (pt)/pt (step 2). 

The dashed curves show the result of steps 3 and 4 — 
modeling the data with boosted soft component S' NN 
(Levy distribution). The solid curves show the RHS of 
Eq. (f2"5|) with the first bracket replaced by 2 /n part to form 
7t(l - /3*)/2T 2 • f(y t ) 2/n par t Ay t2 /02 (?/*)• The minimum- 
bias data used in this analysis correspond to mean par- 
ticipant pathlength v ~ 3.5 as noted in the figure. In the 
product pox{yt) ■ V2Y(yt) of Fig. [H] the correspondence 
proton « Lambda is made for V2(pt) data to estimate 
proton quadrupole spectra. 

Subsequent spectrum interpretation invokes the rela- 
tion of the three quadrupole components to correspond- 
ing hard components on y t7r (dash-dotted curves Hnnx 
in Fig. [5]). The hard components for p and K (all hard- 
component modes are at yt{^) ~ 2.7) strongly overlap 
the corresponding quadrupole components, but that for 
pions does not. Such structural details may explain the 
variation of v^ipt) distributions with mass in relation to 
so-called constituent-quark scaling. 
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B. Quadrupole spectrum and soft component 



IX. THE STRUCTURE OF v 2 



Fig. [9] shows data (solid points) from Fig. [8] trans- 
formed to y t {it, K, p) (proper y t for each hadron species) 
with the appropriate Jacobians. The common left edge 
reveals monopole boost Ay m — 0.6. From Eq. (|25|) the 
form of the data is oc p'Jpt ■ f(yt) • Pi{yt\ Ay t0 ), the last 
factor being the quadrupole spectrum. 
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FIG. 9: Spectra from Fig. |8]transformed to proper yt for each 
hadron species. The dotted curves are soft components from 
respective single-particle spectra for comparison. The promi- 
nent feature is the common edge at yt ~ 0.6, implying that 
the three hadron species originate from a common boosted 
source. The hadron abundances and spectrum shapes are the 
same as the single-particle spectrum soft components. 



The quadrupole spectrum for each hadron species can 
be modeled with the same form of Levy distribution used 
for the soft component of the single-particle spectrum. 
Also plotted in Fig. [5] are soft components SNNx{yt) 
from the single-particle spectra for three hadron species 
(dotted curves). The dashed curves through data points 
are A/T 2 -p'Jptlt^ - Pt) ■ S' NN {y t - Ay m ;T 2 ,n 2 ), with 
factor A and monopole boost Ay t Q common to the three 
species. Levy T 2 and n 2 parameters have been opti- 
mized for each quadrupole spectrum. The factors are 
A/T 2 ~ 0.005/(0.1 GeV) - 1/20 GeV" 1 . The descrip- 
tion of data is good. The solid curves are the same but 
with factor p' t /pt7t{^ — fit) (Fig-Hi right panel) removed, 
revealing the undistorted shapes of p%[yu Ayto). Com- 
parison with the single-particle spectra (dotted curves) 
reveals the similarities of the single-particle soft and 
quadrupole hadron sources. 



In Sec. IVIII v 2 was factored, and in Sec. IVIIII v 2 was 

represented by the combination of a boosted soft compo- 
nent S' NN and two-component single-particle spectrum 
2/n pQr t • po- The full expression with V2/2?r = pa v 2 is 



2 V 2 (y t 



P't 



A 



2tt 



7t(i-A)r 2 
_Jt 

.7t(l-A) 
f{yf,Ay t0 ,Ay t2 ) 



S' NN (y t -Ay t0 ;T 2 ) (26) 
Ay t2 p 2 (y t ;Ay t0 ), 



''part 



where 2/n part ■ p 2 has been modeled by AS' NN {yt — 
Ayto',T 2 ,n 2 ), a boosted soft component (Levy distribu- 
tion) with reduced amplitude and modified shape param- 
eters (T, n). 

The structure of v 2 (p t ) can be understood entirely 
in terms of two yt-dependent factors — p' t and the spec- 
trum ratio. The two-component spectrum model is 
2/npart ■ Po = S NN (y t ;T ) + vH AA (y t , v). The shape of 
the single-particle spectrum in the v 2 denominator varies 
strongly with centrality (v) due to evolution of the hard 
component [lH ]. 



A. p' t = pt in the boost frame 

Fig. [10] (left panel) shows v 2 (pt) data with mass order- 
ing at small p t attributed to hydrodynamic flow 0, H, 0] . 
The mass dependence of v 2 (pt) at small pt is determined 
entirely by p t in a frame boosted on y t by Ay t . p t in the 
boost frame is defined in the lab frame by 



p' t = m sinh(y t - Ay t ) = j t (p t - (3 t m t ) 
= m t 7 t {tanh(y t ) - tanh(Ay t )}. 



(27) 



p' t vs pt in the lab frame is shown in Fig. [TO] (right 
panel). As noted, the curve intercepts are located at 
pto = mo sinh(Ayt). Comparing the two panels it is 
apparent that the "mass ordering" of v 2 attributed to 
hydrodynamics is produced entirely by the kinematic re- 
lation between p' t and pt determined by monopole boost 
Ay t Q. The mass ordering alone does not determine what 
physical mechanism caused the boost. 



B. Spectrum ratios 

Variation of v 2 (p t ) relative to the p' t trend (especially 
above 0.5 GeV/c) is determined by ratios of quadrupole 
to single-particle spectra. Fig. [TT1 shows spectrum ratios 



S NN(Vt - Ay m ;T 2 ) ^ 
S NN (y t ;T ) + v H AA (y tl is) p 



(28) 
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FIG. 10: Left panel: V2{pt) data from Fig. \T\ (left panel) are 
repeated for comparison. The curves through data are from 
the present analysis. Right panel: Relation between pt in the 
boost frame and lab frame for three hadron masses. Intercepts 
on the abscissa are at p t o values defined in the text. 



in two plotting formats for three hadron species (pions 
- solid, kaons - dashed, Lambdas - dash-dot). The nu- 
merator of Eq. (f2"g| appears in Fig. [5] (plotted on pion 
yt) as solid curves in the form A/T 2 ■ S' NN (y t - Ay t0 ; T 2 ) 
[oc quadrupole spectra P2(yt)]- The denominator appears 
in that figure as dotted curves [oc two-component single- 
particle spectra poiiJt)}- Those spectrum models describe 
V2(pt) and spectrum data within published errors. 




p t (lab) (GeV/c) 



2 3 4 
y,(7i,K,p) 



FIG. 11: Ratio of quadrupole (boosted soft component) spec- 
tra to single-particle spectra plotted on pt (left panel) and 
proper yt for each hadron species (right panel; solid, dashed, 
dash-dot) for minimum-bias Au-Au collisions (y ~ 3.5). Spec- 
trum ratios with no hard component (dotted, v = 0) are rel- 
evant to a hydro description. Ratios for no monopole boost 
(dashed, Ayto = 0) are also plotted. The dominant role of the 
hard component in V2(pt) is most evident in the right panel. 

In Fig. QT] the dashed curves at lower left represent 
spectrum ratios for no monopole boost (Ay t o = 0) and no 
hard component (y = 0), which then directly relate the 
quadrupole boosted S' NN shape to the soft-component 
shape of the single-particle spectrum. The ratios are de- 
fined to obtain unit magnitude at yt = Pt = for those 
conditions. Quadrupole spectra inferred from V2{pt) data 
are significantly narrower for each hadron species, i.e. 
substantially "cooler" than the single-particle soft com- 
ponents (T2 ~ 0.7 To). The dotted curves at upper left in 
each panel result from "turning on" the quadrupole boost 
Ay*2 = 0.6 but not the hard component (minijets) in the 



denominator of Eq. (|28p iy — * 0) to form soft-reference 
spectrum ratios. 

The v = 3.5 spectrum ratios (solid, dashed, dash-dot 
curves) include the hard component in the denomina- 
tor. Deviations from the v — (dotted) soft-ratio ref- 
erence curves (arrows) demonstrate the dominant role of 
the hard component in distorting spectrum ratios, and 
therefore V2(pt), over most of the Pt/yt range. The dis- 
tortion is not a consequence of "nonflow" as defined in 
conventional flow analysis (discussed in Sec. IXIV[) . It is 
inherent in the V2{pt) definition as a ratio and results 
from the minijet contribution to its denominator. 

The mass systematics for the quadrupole and soft com- 
ponents are simplest to describe in Fig. QT] (right panel). 
From Eq. (f26|) spectrum ratios are oc V2{pt)/p't- Fig. [Til 
(right panel) is therefore comparable to V2/pt data in 
Fig. Q] (right panel) . Soft-component or "thermal" spec- 
tra (y = 0) on specific hadron y t vary approximately as 
exp{— TOo[cosh(y t ) — 1]/T}. The spectrum widths then 
vary as 1/mo, as illustrated by the dashed curves at 
lower left for Ay t Q = 0. For nonzero Ay t o (~ 0.6) the 
same trend holds for increases above unity of the soft- 
component ratios [y = 0), as illustrated by the dotted 
curves. Deviations of data ratios from the dotted curves 
(arrows) are determined by the relation between soft and 
hard spectrum components for different hadron species. 



C. Role of the hard component in V2(pt) 

Because the modes of the hard components for three 
hadron species are located near 1 GeV/c (reflecting a 
common underlying parton spectrum) spectrum ratios 
for more-massive hadrons are more affected at smaller 
yt- The quadrupole and minijet peaks for kaons and 
Lambdas/protons occur at nearly the same position on 
Pt, whereas the peaks for pions are significantly sepa- 
rated, as in Fig. [S] Thus, the pion ratio is affected only 
above about 0.5 GeV/c, but the kaon and Lambda (pro- 
ton) spectrum ratios are dominated by the hard compo- 
nent over all pt /yt , explaining the mass dependence of the 
sequence of downturns (arrows) from the dotted curves 
in Fig. QT] (right panel). 

In effect, the hard component interacts through ra- 
tio i>2 with the quadrupole component and soft spec- 
trum component. There is a complex numerical inter- 
play between hadrons from a boosted source, longitudinal 
participant-nucleon fragmentation and transverse parton 
fragmentation. The correspondence below y t ~ 1.5 arises 
because the boosted hadron source apparently produces 
hadron species in the same abundance ratios as the N-N 
soft component (cf. statistical model [23, HH). The main 
difference is a smaller slope parameter T 2 ~ 0.7 To. This 
exercise demonstrates the importance of interpretable 
correlation measures and plotting formats. Quadrupole 
and minijet systematics should be studied within inde- 
pendent analysis contexts. In «2(pt) they arc maximally 
confused. 
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X. QUANTITATIVE COMPARISONS 

The description of quadrupole spectra in Sec. |IX] can 
be related quantitatively to other spectrum features and 
trends. The goal should be the simplest and most com- 
prehensive description of all single-particle and correla- 
tion structures. 



at 200 GeV we have 
2 Ay t2 n ch2 



tpart 



2tt 



2 Ay t2 n ch2 

T-part n NN 



2T- 

^0mAhvn NN /2T: (33) 

Pt 
0.028 
0.016, 



A. Relation to the soft-component spectrum 

From Fig. [8] we can infer that 

2 v2x(Pt) „ 0-005 , yoqi 
Pox — ~ -7^b NNX (y t -Ay t0 ;T 2 ,n 2 l29) 



n P art Pt T 2 

for each hadron species X. From Eq. (|25|) we then have 

2 * 2-0.005 . . Wo „, 
Ay t2 p 2 x ^ — -tt S NNX {yt ~ Ayto) (30) 

"part 7t(l-ft) 

~ 0.016 S' NNX (yt-Ayto), 

since 7t(l — /3t) ~ 0.6 and T 2 ~ 0.1 GeV for minimum- 
bias data. Thus, N-N soft particle production from a 
common boosted source describes the quadrupole in A- 
A minimum-bias collisions. 



B. Relation to soft and hard hadron yields 



where the last line applies to the minimum-bias case, with 
v ~ 3.5 and e ~ 0.3. From Eq. (f30|) with n so ft ~ unn 



2 Ay t2 n 



ch2 



Ipart "AW 



0.016 



(34) 



for minimum-bias collisions, which is consistent. An in- 
dependent method is required to place limits on the ab- 
solute multiplicity n c h2 of the quadrupole (cf. Sec. lXIIl| ). 
For comparison, the full-spectrum and hard-component 
integrals relative to the N-N multiplicity are 



rich 



Tlpart "AW 
2 "/lard 
"part "AW 



1 + 0.1(^-1) 

O.lZA 



(35) 



Roughly, 2/n part n haTd ~ 0.25v. 

The factors in the second line of Eq. (|33| can be in- 
terpreted in terms of quadrupole emission as follows: 
The final-state quadrupole moment (hadron pair yield) 
oc (Ay t2 n c h 2 ) 2 goes as interaction length v times the 
initial-state quadrupole moment oc (n part /2 ■ e op t) 2 . 



Minimum-bias spectrum data do not provide informa- 
tion about centrality dependence. In [15| the centrality 
trend for p t -integrated v 2 was inferred from v 2 {4} data 



1 



Po (2^) 2 
2 V 2 
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= 0.0045 e 
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(31) 



= . / 0.0045 1/- 



2 n ch 

ripart 27T 



The per-participant hadron quadrupole density squared 
oc v, whereas for minijets the per-participant hadron 
fragment density oc v, possibly revealing the difference 
between quadrupole radiation and parton scattering [l5| . 

Integrating Eq. ([23)) . with n c h2 the integral of 
quadrupole spectrum p 2 over p t and one unit of rapidity, 
gives 



V 9 



p' t 1 2 Ay t2 n ch2 



n part 2ir e 2T 2 e n par t 



2n 



(32) 



where Ay t2 n c h 2 is the effective number of quadrupole 
hadrons in one unit of rapidity. Given the v 2 (pt) data 
n ch2 could be a large number with weak boost modu- 
lation Ay t2 or a small number with strong modulation. 
Invoking 2/n part -n ch ~ rijviv(l + 0.1[i' — 1]) ~ n NN = 2.5 



C. Relation to minijet pair correlations 

The hard-component spectrum yield should relate to 
observed changes in pair correlations associated with 
minijets. The factor 



n p art 



1 dn 

2tt dn ^° 



0.4(1 + 0.1[u - 1])(36) 



times Ap/y/p Te f (per particle) gives the number of cor- 
related hadron pairs per participant nucleon pair. The 
density of correlated fragment pairs per participant pair 
is therefore 



n p art 



-Ap(?7A,<M) = 0.4(1 + 0.1 [i/ - 1]) 



Ap 



(37) 



The integral of Ap^/p re f over the same-side minijet peak 
increases by 6x over binary-collision scaling. The in- 
crease in number of correlated pairs per participant over 
binary collision scaling for central collisions is therefore 
1.5-6 = 9x . That means the number of hadron fragments 
per parton pair increases by 3x from ~ 2.5 to ~ 7.5. 

For binary collision scaling of N-N collisions the frac- 
tional increase of hadrons which are correlated minijet 
fragments from peripheral to central Au-Au should be 
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8%. It is observed to be 25% due to the factor 3x derived 
above. The fraction of total hadrons which is correlated 
fragments is thus 25/125 = 20%. 

The hard-component fraction of single-particle spectra 
is ~ 30% for central Au-Au collisions. The hard compo- 
nent of spectra appears entirely in narrow structures for 
pions and protons corresponding to the same boost [l4[ • 
The 13% difference is therefore intimately connected with 
the anomalous boost component, which is in turn con- 
nected with minijet broadening. 
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XI. QUADRUPOLE VS HYDRO THEORY 

Using this detailed description of V2(j>t) structure we 
can explore the relationship of the quadrupole compo- 
nent to theory predictions. Hydrodynamics is tested by 
the quadrupole component in three ways: the source 
boost distribution, the apparent "temperature" T% of 
the quadrupole source and the abundances of the sev- 
eral hadron species produced by the source. Parameters 
of the spectrum soft component, such as temperature To, 
are obtained from measured single-particle spectra. 

In Fig. [1] dotted curves in the two panels represent 
viscous hydro calculations for pions [l], H3] which have 
been interpreted to support claims for formation of a 
low-viscosity medium, possibly a "perfect liquid." The 
theory curves represent the limiting case of zero viscosity. 
Although the curves in the left panel appear to describe 
the pion data, transformation to the right panel reveals 
that theory and data are actually substantially different. 



FIG. 12: Left panel: Ratio of quadrupole component (boosted 
soft component) to spectrum soft component (both Levy dis- 
tributions) for three hadron species illustrating typical hydro 
mass dependence of the ratio. The dotted curve includes fac- 
tor Pt/p*7t(l — Pt) for comparison with the right panel. Right 
panel: v^jpt data compared to hydro theory curves as de- 
scribed in the text. Different boost distributions are apparent 
at lower left. 



mass sequence on proper y t . The ratios correspond to 
the dotted curves in Fig. [TT] (right panel), but the am- 
plitudes are adjusted so all curves start at the same ini- 
tial value. The dotted curve, including additional factor 
p' r /"ftO- + Pt)pt, is within a constant factor the solid curve 
in the right panel. 



B. Data comparisons with hydro theory 



A. Spectrum ratios 

Typical hydro calculations do not include a hard com- 
ponent in the denominator of 1*2 (p*)- Thus, the spectrum 
ratio in Eq. (|28p becomes 

S' N N(yt - Aj/to) ^ exp{-7 t (m t - /3 t Pt)/T 2 } 
SNN(yt) exp{-m t /T } 

~ exp{[l/T - 7i (l-A)/T 2 ]m t }, 

where the boosted and single-particle soft components 
are approximated by exponentials with possibly differ- 
ent temperatures, and 74, fit are determined by monopole 
boost Ay to- In a quantitative analysis the soft compo- 
nents are modeled by Levy distributions with possibly 
different indices n. The ratio is ~ an exponential with 
positive constant determined by the temperature differ- 
ence and the monopole boost, as in the second line. 

In Fig. [T2] (left panel) the solid, dashed and dash- 
dot curves show spectrum ratios as in Eq. I|38p for three 
hadron species using spectrum soft components (denom- 
inators) which properly describe single-particle spectra 
(i.e., Levy distributions) and boosted components (nu- 
merators) with the same temperatures for simplicity and 
boost Ay to = 0.6 as for V2(pt) data, showing the expected 



In Fig. [T2] (right panel) the data from Fig. Q] (right 
panel) are repeated for comparison (the three curves 
through the data points are from this analysis). The 
lower dotted curve is the zero-viscosity hydro prediction 
from [l[ . The solid curve following the hydro dotted curve 
at larger y t is Bp'Jpt ■ S' NN {y t - Ay m ) / S NN (y t ) (oc the 
dotted curve in the left panel). The soft-component ratio 
S' NN /S NN is defined in Eq. (3SJ|, with T = 0.14 GeV, 
T 2 = 0.095 GeV and Ay m = 0.6. Factor p'Jp t has been 
added to incorporate the form plotted in Fig. [7] (right 
panel) appropriate for the viij>t)lpt ratio. B is adjusted 
to match the hydro (lower dotted) curve at larger y t . 
Agreement of the shapes is good except near the origin 
where the boost distributions differ. 

The dashed and upper dotted curves in the right panel 
are 2.7x the solid and hydro curves. The dashed curve 
describes the data for three masses well in the smaller-y t 
region where the hard component does not dominate the 
variation, as expected from Fig. [TT] (right panel). The 
zero- viscosity hydro curve for pions from [l[ thus under- 
predicts the V2 magnitude by 2.7x in p t < 0.5 GeV/c 
(y t7T < 2). This exercise is intended to demonstrate how 
hydro theory can be better tested in the plotting format 
of Fig. [T21 in particular the predicted boost distribution. 
The validity of specific hydro theories is the subject of 
subsequent analysis. 
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C. Comparison of boost distributions 

In Fig. [TJ] (right panel) the structure in y t < 1.5 is 
possibly the first direct comparison of boost distribu- 
tions from data and hydro theory. Boost comparisons 
provide essential tests of hydro and the expanding bulk 
medium scenario for heavy ion collisions. Boost details 
are strongly suppressed in plots of V2{pt) vspt- The data 
for kaons and Lambdas in the form V2(pt)/pt on proper 
yt clearly contradict the hydro prediction below y t ~ 1.5 
where the boost distribution dominates. 

The data indicate a narrow boost distribution centered 
at Ayto = 0.6. The hydro prediction suggests a broad 
distribution starting at yt = and roughly consistent 
with Hubble expansion of a bulk medium. An essential 
requirement for any theory of the quadrupole component 
is an explicit boost distribution compared with accurate 
data. Inferences of small (or any) viscosity from compar- 
isons of V2(pt) data with the lower dotted curve [l| are 
not justified. 

In Fig. [T2] (right panel) there is the suggestion of real 
"scaling" derived from a common boosted hadron source, 
the correspondence of the three hadron species in this 
plotting format below yt ~ 1.5. Above that point the 
curves deviate from the hydro hypothesis by large fac- 
tors determined by spectrum hard components (parton 
scattering and fragmentation) . Each hadron species devi- 
ates from the universal hydro curve at a point depending 
on its mass, revealing interaction of soft components with 
the universality of the underlying parton spectrum onp t . 

XII. QUADRUPOLE MODEL UNCERTAINTIES 

The quadrupole spectra in Figs. [5] and O the main re- 
sults of this analysis, were obtained as simple combina- 
tions of previously measured data. As such, uncertainties 
indicated by error bars are propagated from the original 
published errors, but spectrum parameters inferred from 
those data possess unique uncertainties to be estimated. 

The common left edges in Fig. [9] taken together de- 
termine Ay to = 0.6 ± 0.05. The boost distribution ap- 
pears to be narrow (r.m.s < 0.1) even though this is a 
minimum-bias centrality sample, but the data are too 
sparse in that region to provide a better width estimate. 

The spectrum shapes near the left edges determine 
T 2 ~ 0.1 ± 0.005 GeV (~ 0.09 GeV for pions and - 0.11 
GeV for protons). The shapes further out on the tails of 
the distributions determine Levy exponents n 2 ~ 15, but 
the shapes are also influenced by /(y t ; Ay t0 , Ay t2 ) > 1. 
Since the spectra in Fig. [5] are uncorrected for that 0(1) 
factor the T 2 estimates should be taken as upper limits 
and the Levy index n 2 estimates as lower limits. 

The largest uncertainties apply to the estimates of ab- 
solute quadrupole yields. In Fig. [9] the quadrupole am- 
plitudes at spectrum left edges determine the relative to- 
tal yields, which correspond well to the single-particle 
spectrum soft-component relative yields (dotted curves 



in that figure). Yield uncertainties from /(yt), which 
mainly affects the spectrum tails at larger y t , are small. 

The quadrupole spectra contain a common factor Ay t2 
which is the major source of uncertainty in estimating 
the total spectrum yields, as discussed in Sec. IXIII1 The 
uncertainty in the absolute yields is less than a factor 2 x , 
sufficient to determine that the quadrupole component is 
at most a small fraction of the total particle yield. 

Quadrupole y t spectra observed directly are simple, de- 
scribed by a few parameters and very similar in shape to 
single-particle spectrum soft components, albeit boosted 
on y t . However, when coupled to two-component spectra 
via ratio ^(pt) "elliptic flow" the data become arbitrarily 
complex and essentially uninterpretable. 



XIII. QUADRUPOLE ABSOLUTE YIELDS 

The quadrupole absolute yield n c h2 (in one unit of 
rapidity) can provide definitive model tests. But from 
V2{pt) data alone there remains an ambiguity in the prod- 
uct Ay t 2 n c h2 ■ The ambiguity is reduced by the edge 
of the quadrupole spectrum near y t ~ Ayto, evident in 
Fig. UJ (right panels), but accurate data in that region are 
difficult to obtain. 

The single-particle spectrum hard component, de- 
scribed as minimum-bias parton fragmentation to mini- 
jets, suggests a solution. The hard component was iso- 
lated by a combination of techniques: correlation anal- 
ysis of several hadron charge- sign combinations [HI, H|| , 
n c h dependence of p-p spectra [28j , v dependence of Au- 
Au spectra [141 ] and comparisons with the systematics of 
fragmentation functions from e + -e~ collisions [37l ]. 

In this section comparisons of spectrum shapes and 
centrality trends for integrated yields are combined to es- 
timate the absolute quadrupole yield. Quadrupole spec- 
tra extracted from vi(j>t) data are compared to single- 
particle Au-Au spectra, and centrality trends of pt- 
integrated quadrupole data are compared with those of 
integrated single-particle spectrum structures. 



A. Lower limits from the soft-component model 

The boosted soft-component model of quadrupole data 
(dashed curves) in Fig. [5] can be used to provide a lower 
limit to quadrupole spectra. The dashed model curves 
are defined according to Eq. (|26|) by 

— ^ = { * a , \ p' NN (y t ~ A ytoX 39) 

npart Pt lPtlt(l ~ Pt) J T 2 

with A ~ 0.005 and A/T 2 ~ 0.05. From Eq. (|24|) 

Po«2 = ^^ - ^rf(yt)Ay t2 p 2 (y t ;Ay t0 ). (40) 
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A requirement of positive-definite boosts implies Ay t 2 < 
Ayto ~ 0.6 and we obtain 



<>part 



-Pi 



2 AS', 



NN 



lt (l-P t )Ay t2 



> 0.025S" A 



NN- 



(41) 



The lower limits are represented by the solid curves in 
Fig. [T3J which are about 0.5 x the solid curves in Fig. [5] 
because 0.025T 2 /A - 0.5. 




FIG. 13: Lower (solid) and upper (dotted) limits on 
quadrupole spectra, the latter obtained by direct comparison 
to single-particle spectra (dash-dot curves). 



B. Upper limits from single-particle spectra 

Loose upper limits on pi(yt-K\ Ay to) can be estimated 
by direct comparison with the full single-particle spec- 
tra. Factor f(yt) is ignored, since only the most promi- 
nent spectrum aspects at smaller y t matter, especially the 
edges of the boosted distributions. In Fig. [T3] the rough 
upper limits on quadrupole spectra (dotted curves) are 
determined by the condition that they not exceed 10% 
of the measure single-particle spectra at any point. The 
upper limits are then only a factor 2x the lower limits. 



C. Upper limits from spectrum residuals 

Tighter constraints can be established by comparing 
quadrupole spectrum shapes to residuals of a compar- 
ison between single-particle spectrum data and a two- 
component spectrum model [14]]. Fig. [TJ] (left panel) 
shows minimum-bias quadrupole spectra for pions and 



protons (thick solid curves) compared to the residu- 
als from minimum-bias (y ~ 3.5) single-particle pion 
and proton spectra (thin curves and open symbols) ob- 
tained by comparing spectrum data with an accurate 
two-component (soft-plus-hard) model [3]. The pro- 
ton residuals peak corresponds to the p/7r ratio "puz- 
zle" [2l], I23I ? ] . The even- larger pion residuals peak was 
not previously noted. 



0.008 




FIG. 14: Left panel: Minimum values of quadrupole spectra 
compared to minimum-bias spectrum residuals relative to the 
two-component model [3]. Right panel: Expected central- 
ity dependence of the quadrupole integral 2/n par t Ay t 2n c h2 
compared to that for the spectrum hard component. 

The quadrupole (thick solid) curves have the minimum 
amplitudes determined above by assuming a positive- 
definite radial boost Ay t 2 ~ Ay to. The spectrum residu- 
als appear inconsistent with any increase in quadrupole 
amplitudes beyond the minimum. 



D. Upper limits from centrality dependence 

To confirm the upper limits from minimum-bias spec- 
trum residuals detailed centrality dependence of spec- 
trum structure can be compared on proper hadron ra- 
pidity for each species. In Fig. [14] (right panel) the 
solid curve shows the quadrupole centrality dependence 
in Eq. g2) 



Ay t 2 n, 



ch2 



n NN 0.028eVv 



(42) 



from the analysis in [l5[ ■ The centrality dependence of 
the hard component is also sketched for contrast. The 
minimum-bias v^ipt) data used in this analysis corre- 
spond to the maximum of product Ay t 2 n c h2 on central- 
ity. The product should decrease strongly for more cen- 
tral collisions compared to hard-component structure. 

In Fig.[l5]the lower limits of minimum-bias quadrupole 
spectra for pions and Lambdas (thick solid curves) are 
compared to the centrality variation of single-particle 
spectrum residuals on proper hadron y t . The common 
monopole boost Ay t o = 0.6 for the pion and Lambda 
quadrupole curves is apparent. 

Pion and proton spectrum residuals for five centrali- 
ties are indicated by the thin curves of different types 
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FIG. 15: Lower limits on quadrupole spectra in the form 
"2/n-part p2(yt) (solid points and curves) compared to single- 
particle spectrum residuals (open points and thin curves of 
various styles) for five Au-Au centralities relative to a two- 
component reference [lj. Soft components SNN{yt) provide 
a reference. The dashed curves corresponding to hadrons from 
a boosted source with Ay t ~ 1.1 and substantially smaller 
slope parameters T than the quadrupole components suggest 
a possible mechanism for the residual spectrum structure. 



and open symbols. Centralities correspond to the points 
in Fig. [TJ] (right panel). Comparison with the boost sys- 
tematics of the quadrupole spectra suggests that the pro- 
ton and pion spectrum residuals may also arise from a 
common boosted source, but with quite different boost 
distribution. The thick dashed curves sketch a com- 
mon boosted source with Ayt ~ 1.1 for the two hadron 
species. In that context the boost distributions appear to 
be strongly centrality dependent. The peak modes move 
to smaller rapidities for more central collisions. 

The solid points are quadrupole data from Fig. [9] with 
minimum amplitudes as in Fig. [H] (left panel), kine- 
matic factor Pt/ptTtO- ~ fit) removed and multiplied by 
0.025 T 2 /A ~ 0.5 according to Eqs. (3SJ) and ©J. The 
arrow at left indicates a Lambda v-i datum which is neg- 
ative but consistent with zero. The agreement between 
points and solid curves demonstrates that boosted Levy 
distributions describe the quadrupole spectra well. 

Comparisons of quadrupole data and spectrum resid- 
uals indicate that increasing the quadrupole magnitude 
beyond the lower limit would strongly conflict with the 
residuals. The structure and centrality dependence of the 
proton residuals appear to be consistent with the extrap- 
olated centrality evolution of the lower-limit quadrupole 
component. The centrality dependence of the structure 
at the left side of the proton residuals peak may arise 



from the quadrupole contribution. The pion compari- 
son is indeterminant because of a lack of spectrum data 
below y t ~ 1.7. This detailed comparison is limited by 
sparse data and data uncertainties but suggests that the 
upper limit on the quadrupole component is consistent 
with the lower limit. A more precise statement requires 
improved spectrum and quadrupole data at small y t . 

E. Discussion 

The small upper limit (~ 5%) on the fraction of final- 
state particles participating in the azimuth quadrupole is 
certainly counter-intuitive. The conventional scenario for 
more-central RHIC collisions is that almost all hadrons 
emerge from a common partially-thermalized medium 
supporting radial and elliptic flow. The definition of 
V2(pt) implicitly relies on the assumption that the single- 
particle spectrum in the denominator is the same as the 
quadrupole spectrum contained in the numerator com- 
bined with other factors. 

A central message of the present analysis is that the 
quadrupole and single-particle spectra are not the same, 
that the former is boosted significantly and the boost is 
not shared by the single-particle spectrum (radial flow 
is negligible [14[). Distinctions between spectrum shapes 
are the basis for estimating what fraction of the final 
state actually carries the azimuth quadrupole structure. 

Since i>2 measures the product of the true momen- 
tum asymmetry and the fraction of particles carrying the 
quadrupole, small values of V2 plus the conventional as- 
sumption about a flowing bulk medium suggest that the 
momentum asymmetry is rather small (few percent), and 
therefore can be explained by hydrodynamic response to 
initial pressure gradients. 

Comparing the reconstructed quadrupole spectrum 
with the single-particle spectrum reveals that the 
quadrupole fraction is actually small, and therefore that 
the momentum eccentricity for that small fraction of par- 
ticles is large (near the upper limit defined by a require- 
ment of positive-definite boost). That conclusion is not 
inconsistent with any previous i>2 measurements, only 
with a priori expectations within the hydro context. 



XIV. QUADRUPOLE VS NONFLOW 

Fig. [IT] demonstrates that the hard component of the 
single-particle spectrum present in the denominator of 
i>2 severely distorts the structure of v^ipt) above about 
0.5 GeV/c for pions, kaons and protons. The hard- 
component distortion should be distinguished from possi- 
ble nonflow distortions also due to minijets but appearing 
in the numerator of v 2 . 

Nonflow is dominated by minijet angular correlations 
misinterpreted as azimuth quadrupole correlations by 
conventional ID flow analysis methods [l5j. Minijet cor- 
relations and the spectrum hard component have a com- 
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mon source: hadron fragments from a minimum-bias 
scattered parton spectrum [la ]. The combination pro- 
duces large uncertainties in the interpretation of V2 data 
above 0.5 GeV/c. In this section I consider the mecha- 
nism and consequences of nonflow contributions to «2 (j>t) ■ 



A. Nonflow and the hard spectrum component 

The structure of V2{2}(p t ) (ID azimuth correlations) 
including nonflow is described schematically by 



v 2 {2}(pt,is) oc p' t 



i AS NN(Vt- Ay ) + nonflow 



SNN{yt) + vH AA {y tl v) 



.(43) 



"Nonflow" is dominated by the m = 2 azimuth Fourier 
amplitude of the same-side mini jet peak (jet cone) in an- 
gular correlations on (rj,(f) [IH, [38|. The relative mag- 
nitude of the nonflow term in vi depends in part on the 
analysis method and spectrum structure. At larger p t for 
pions and all pt for less-abundant hadrons the V2{EP} 
(event-plane) method is typically employed to accommo- 
date smaller particle yields. V2{EP} ~ ^2{2} is maxi- 
mally sensitive to minijets [Til ]. In contrast, 2D angular 
autocorrelations on (77, <j)) can be used to separate minijet 
and quadrupole components accurately [HI, [H| . 

Hard component vHaa (yt , v) is the angle- integrated 
minijet fragment spectrum, whereas nonflow is a Fourier 
component of the minijet same-side peak on azimuth. 
Thus, minijet contributions in numerator (nonflow) and 
denominator (hard component) of 1^2 (p*) are directly re- 
lated. However, the nonflow contribution has its own 
substantial pt dependence relative to the spectrum hard 
component (i.e., minijet yield) because the Fourier ampli- 
tude of the same-side peak depends on the peak shape (77 
and <f> widths), which varies strongly with parton energy 
scale (as determined by the selected hadron fragment pt). 



B. V2(pt) trends at larger pt 

In Fig. [16] quadrupole spectrum components inferred 
from this analysis are compared to soft and hard single- 
particle spectrum components for pions and protons from 
minimum-bias Au-Au collisions. With increasing p t there 
is competition between the tails of the quadrupole spec- 
tra and hard-component spectra. The latter completely 
determine V2(pt) trends at larger p t . 

Above pt ~ 2 GcV/c (yt n ~ 3.3) the soft-component 
and quadrupole spectra are dominated by the hard 
component (parton fragments) 0, H3- By definition 
V2{i\{pt) represents the m = 2 Fourier component of 
all azimuth correlation structure. Thus, at larger p t 
V2{i\{pt) is simply the ratio of the Fourier amplitude 
of the same-side minijet peak (jet cone) to the spectrum 
hard component (minijets). v^ipt) then follows a nearly 
constant trend on p t described as "saturation." Varia- 
tion of the ratio (modulo prefactor p' t ) is dominated by 




FIG. 16: Comparison of the quadrupole component with soft 
and hard spectrum components for pions (left panel) and pro- 
tons (right panel). 



changes in the minijet peak shape with parton energy 
scale. There is no required relation to a reaction plane. 

The centrality variation of V2{pt) at larger p t should 
be dominated by parton energy-loss effects, including r\ 
broadening of the same-side peak in the numerator (non- 
flow) [38] and reduction of the hard-component tail in 
the denominator ("jet quenching") [TJ]. Thus, in a con- 
ventional flow analysis the region above 0.5 - 1 GeV/c 
is already strongly distorted by the spectrum hard com- 
ponent and thus difficult to interpret. The region above 
3-4 GeV/c provides no information about the quadrupole 
component, whether hydro or simple boost phenomenon. 



XV. QUADRUPOLE VS SCALING 

Certain scaling relations are inferred for V2(pt) to 
support claims that "elliptic flow" is a hydrodynamic 
phenomenon manifested by a thermalized bulk partonic 
medium [3, H, that hadronization from the medium 
proceeds via coalescence/recombination of constituent 
quarks [Til . [3^ |. relating to a similar model of certain 
spectrum features (e.g ., the anomalous p/ir ratio at in- 
termediate Pt)(2ll. |22|. [23l| . The overarching conclusion 
from V2{pt) scaling is that sQGP (a thermalized, small- 
viscosity bulk partonic medium) has been formed. In 
this section claims of constituent-quark and other forms 
of v^ipt) scaling are re-examined in the context of the 
present analysis. 



A. V2 scaling observations 

Arguments in favor of a locally-thermalized pre- 
hadronic bulk medium evolving according to near-ideal 
hydrodynamics include: 1) the minimum-bias multiplic- 
ity distribution form is independent of system size, 2) 
hadron species abundances follow a statistical model, 3) 
large i>2 values reveal rapid thermalization and early pres- 
sure gradients common to all hadron species and incom- 
patible with hydro evolution of hadrons (e.g., D and 
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(f) meson V2 data are interpreted to imply thermaliza- 

tion) m. 

Scaling relations invoked in flow studies are inter- 
preted to buttress the above arguments. Scaling rela- 
tions involve combinations of V2(j>t) data, m t (transverse 
mass) and n q (constituent quark number). The mass 
dependence of V2(pt) at small pt is attributed to hy- 
drodynamics. Constituent quark scaling expressed by 
v i{Pt) = nsL v 2( n qPt)i with n q — 2 for mesons and 3 for 
baryons [ll| is interpreted to imply hadronization from 
a thermalized partonic medium. 

In [39j | a universal scaling of the combination 
v 2(Pt) I e n q vs "kinetic energy" (m t — mo)/n q was claimed 
over a broad range of centralities, strongly suggesting 
formation of a thermalized partonic medium. However, 
other measurements disagree with the claimed universal 
centrality trend [26j |. Universal scaling results are also 
claimed for 30-70% centrality, but one can ask when is 
the system not in equilibrium? For what circumstances 
do such scaling trends not hold? What collision systems 
(e.g., N-N) do not thermalize or form a "perfect liquid?" 

The present analysis strongly suggests that most 
hadrons emerge from several nearly-independent QCD 
processes (nucleon or parton scattering and fragmenta- 
tion), but some coupling develops among the processes 
in more-central Au-Au collisions. In Sec. lIXl it was shown 
that there are two shape factors in V2{pt)'- Pt m the boost 
frame and the spectrum ratio S' NN (yt — Ay t o) / po(yt) ■ In 
what follows I consider scaling arguments for each factor. 



B. V2 scaling and p[ 

In Fig. |TD] (right panel) p' t vs p t is plotted. 7 f (1 — 
fit) ~ 0.6 common to three hadron species determines 
all structure. Similar "mass scaling" of V2{pt) is taken to 
imply hydrodynamic flow. But the mass dependence near 
the origin is determined by a single radial (monopole) 
boost Aj/ t o, and there is no indication from such data of 
the actual boost and hadron production mechanisms. 



> 





12 3 4 
2(m t - m )/n q (GeV/c) 



FIG. 17: Left panel: pt in the boost frame compared to trans- 
verse kinetic energy mt — mo in the lab frame for monopole 
boost value Ayto = 0.6. The curve intercepts are at mto — mo 
as defined in the text. Right panel: The same relations 
rescaled by "constituent quark number" as 2/n q so that me- 
son trends are unchanged. The shift for baryons is indicated. 



In Fig. [17] (left panel) p' t is replotted on m t — m . The 
mass dependence near the origin appears to be reduced, 
but the locations of the curve intercepts are simply given 
by mm " m a = m (cosh[Az/ t0 ] - 1) ~ m (Ay t0 ) 2 /2 on 
m t - ra compared to p t0 = m sinh(Ay t0 ) — m Ay t o 
on pt noted in Fig. 2] (left panel). Consequences of the 
source boost, especially boost distribution details, are 
compressed on m t — too (by a factor 3 for Ay t o ~ 0.6) in 
the pt region most important to the hydro interpretation, 
but the boost is just as accurately determined from the 
data regardless of plotting format. Fig. [1] (right panel) 
clearly provides the best visual access. At larger p t p' t — > 
m% is shifted upward by mo relative to abscissa m t — too. 

In Fig. [T7] (right panel) both axes are scaled by 2/n q 
(factor 2 so axis values remain the same for comparison). 
The consequences are trivial. The intercept at smaller 
nit is reduced by 2/3 for baryons, and the constant ver- 
tical offset mo at larger m t is also reduced by 2/3 for 
baryons. Visual differences between baryons and mesons 
are indeed reduced, but the results are not fundamental 
because the form of V2(pt) at larger p t is not determined 
by hydro or any boost phenomenon (cf. next subsection). 



C. v-2 scaling and spectrum ratios 

Fig- HH (left panel) shows v^{pt) data for three hadron 
species plotted in the same format as Fig. [T7] (left panel). 
This figure can be compared directly with Fig. 4 of [331. 
The dotted curves represent the hydro theory curve [l| 
and 2.7x hydro. Data near the origin follow the p' t sys- 
tematics described above. As in Sec. IIXBI (spectrum ra- 
tios) the turnover of vi(j>t) above 0.5 GeV/c is due to the 
hard component in the V2 denominator. If V2 data do not 
return to zero at larger m t a significant nonflow contri- 
bution is probably present, as discussed in Sec. IXIVl 




m t - m (Ge V/c ) 



2(m t -m )/n (GeV/c*) 



FIG. 18: Left panel: V2(pt) is plotted vs kinetic energy 
mt — mo- The intercept spacing at small p t is reduced by 
a factor 3 x . The lower dotted curve is the zero- viscosity hy- 
dro curve [l[. The upper dotted curve is the hydro curve 
x2.7. Right panel: The left panel except with constituent 
quark scaling in the form 2/n q so that meson trends remain 
unchanged. The shift for baryons is indicated by the arrows. 

In Fig. [TH] (right panel) the n q scaling strategy is used 
to minimize apparent differences between baryon and me- 
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son data, the resulting shifts indicated by the arrows. 
The most dramatic changes occur above 1 GeV/c where 
the data are not relevant to a hydrodynamic mechanism 
or soft processes. In scaling exercises the region above 1 
GeV/c is viewed as dominated by elliptic flow and soft 
hadronization. Scaling trends there are interpreted in 
turn to imply that hadron production is dominated by 
coalescence of constituent quarks. 

The apparent correspondence of data for different 
hadron species left of the vertical dashed lines indicates 
that all-important information about the source boost 
distribution (cf. Fig. [12] - right panel) has been made 
visually inaccessible by a simple transformation. Note 
the limited region of comparison between the upper hy- 
dro (dotted) curve and the data. Generally, comparison 
of different hadron species on p t or m t rather than y t is 
unsuited for hydro (common boost) phenomena. 

V2 data to the right of the vertical dashed lines can 
reveal nothing about hydrodynamic phenomena. The 
v 2{Pt) ratio there is dominated by a complex mixture 
of hard processes (parton scattering and fragmentation) , 
soft spectrum components and quadrupole components, 
with different shape parameters (T, n, etc.) for each com- 
ponent. Sec. IIXB1 reveals that the systematics of Fig. [18] 
are determined by the mass dependence of several spec- 
trum components reflecting soft and hard processes. 

The present analysis demonstrates that quadrupole 
spectra are similar to soft spectrum components (Levy 
distributions) unchanged from N-N collisions. The 
quadrupole hadron production mechanism may well be 
the same as in elementary collisions. Inference of "con- 
stituent quark scaling" from V2 data is prompted by a 
combination of several conventional collision mechanisms 
confused by a poorly-designed correlation measure. 



combined with quark-number scaling (dynamical DoF 
are constituent quarks) "suggest that strongly-coupled 
matter with sub-hadronic degrees of freedom may be cre- 
ated in heavy ion collisions at RHIC" 12j. The present 



analysis is inconsistent with those conclusions. 



B. The fragmentation alternative 

The complementary limiting case is A-A collisions 
modeled by linear superposition of N-N collisions accord- 
ing to the Glauber model, and hadron production by in 
vacuo nucleon and parton fragmentation as in elemen- 
tary hadronic collisions. That linear reference invokes 
two independent fragmentation processes to describe A- 
A collisions: participant-nucleon (soft, longitudinal) frag- 
mentation leading to a soft component of pt spectra and 
minimum-bias large-angle-scattered parton (hard, trans- 
verse) fragmentation leading to a hard component. 

Some aspects of fragmentation can appear thermal, 
even described in part by the statistical model, although 
there is no transport via binary collisions (rescattering) 
in the Boltzmann sense. Fragmentation is a maximum- 
entropy process, the entropy maximization achieved via 
splitting cascades. Parton fragmentation in LEP e + -e~ 
collisions is described to the statistical limits of dat a by 
the beta distribution, a maximum-entropy function [371 ] - 
Deviations from the linear-superposition model in more- 
central A-A collisions could result from a few secondary 
parton interactions. Any non-linearities require careful 
differential study relative to a linear reference (the two- 
component model) [3, [Hj]. The burden should be on 
claims of thermalization to rule out fragmentation as the 
dominant mechanism of A-A collisions at RHIC. 



XVI. DISCUSSION 
A. The conceptual context of elliptic flow 

The conventional elliptic flow context is a limiting case 
in which V2 measurements are interpreted to conclude 
that 1) a monolithic bulk medium produced early in the 
collision is "partonic" (QCD quanta dominate the dy- 
namics), 2) the medium is thermalized rapidly via par- 
tonic rescattering and 3) hadrons emerge late from the 
medium via coalescence of constituent quarks. Multi- 
strange v 2 data for instance exclude slower flow develop- 
ment via hadron thermalization (rescattering) (2o| . 

Spectrum and elliptic flow systematics are used to con- 
clude that "constituent quarks" play a role in hadroniza- 
tion. Recombination (quark coalescence) [TJ HH H2, HII 
reproduces "many features" of hadron spectra in the 
intermediate p t region [1.5,5] GeV/c according to [l2j]. 
Hadronization by quark coalescence is also inferred from 
"scaling" of V2(pt) at intermediate pt, where V2 is said 
to "saturate" at a number apparently oc n q above 1 
GeV/c [l2l ]. Limiting values of V2 (i.e., hydro limits) 



C. Importance of measure design 

Ratio measures applied to RHIC data typically con- 
fuse several collision mechanisms. The Raa spectrum 
ratio mixes soft and hard spectrum components. Most 
of the hard component (at smaller p t ) is obliterated by 
the soft component in a measure intended specifically 
to study parton energy loss [l4|. Ratio t>2(pt) similarly 
mixes soft and hard spectrum components. This analy- 
sis demonstrates that the (soft) quadrupole component 
(boosted source) is severely distorted by the hard com- 
ponent (parton fragmentation) over most of the pt ac- 
ceptance. The benefits of improved measure design are 
suggested by comparison of Fig. Q] (left panel) with Fig. [5] 

Some conventional single-particle spectrum analysis 
also produces misconceptions. The results of monolithic 
"power-law" function [40j fits to multicomponent pt spec- 
tra cannot be interpreted (28|. If the entire spectrum 
below some p t value (e.g., 2 GeV/c) is described by a 
blast-wave model (30j the abundant hard spectrum com- 
ponent (minijets) is injected into the hydro parameteriza- 
tion, confusing parton fragmentation with hydrodynamic 
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(Hubble) expansion [l4| . Better understanding of RHIC 
collisions requires a comprehensive differential approach 
to single-particle and correlation measurements, includ- 
ing comparisons to well-defined references. 



D. Comparison of boost models 

In Sec. IVI CI two radial boost models were described 
and model 2, multiple hadron sources including a boosted 
quadrupole component, was adopted for the present anal- 
ysis. Model 1 is the conventional thermalized partonic 
bulk medium common to all soft hadrons. I reconsider 
the model choice in light of the analysis results. 

Model 1 is essentially the blast-wave model of heavy ion 
collisions [3lj applied by hypothesis to almost all particle 
production. Uniform Hubble expansion is assumed for 
longitudinal and radial boosts. The longitudinal system 
is boost invariant; the transverse boost depends on radius 
(Hubble expansion) and azimuth (elliptic flow). Particle 
emission angle <p p is distinguished from particle source 
azimuth <fi s and the normal to the emission surface (j>b- 
The source pseudorapidity -q (polar angle) is not generally 
the same as the particle longitudinal rapidity y z . 

In the present analysis the soft component, quadrupole 
component and hard component are decoupled. I model 
the quadrupole component by normal emission from a 
cylinder at mid-rapidity and z = 0. The blast- wave 
model simplifies to y z = rj = and <p s — <j) p = </>{,. 
Eq. (11) of [U then becomes the first line of Eq. (Til} 
of this paper. Since this analysis emphasizes qualitative 
study of algebraic structure the simplifications are rea- 
sonable. The boost model of [3l| includes po — > At/to and 
P2 — * Ayi2. However, parameters p relate to a Hubble 
expansion model whereas the Ay t relate to an expand- 
ing cylindrical shell. Ay t Q ~ 0.6 po and similarly for the 
quadrupole. Model 1 is thus a limiting case of model 2. 

The model- 1 expectation is that T and po common to 
most particle production are obtained from p t spectrum 
fits, and p2 is obtained from fits to V2 data. The re- 
sult in [3l| for the monopole component is po ~ 0.9 or 
Ayta ~ 0.55 independent of centrality. Close inspection 
of spectrum fits however reveals that the description of 
the overall spectra over the full pt range is poor, espe- 
cially for pions (the model is fitted to data over a very 
restricted pt interval). Compare with the detailed spec- 
trum description in [lj] in the pt interval [0.2,12] GeV/c 
where no radial boost was required. Model-1 attribution 
of a transverse boost system to the majority of particles 
appears to fail, consistent with the present analysis. 

The approximate quadrupole spectra determined by 
data points in Fig. [5] are obtained from a simple com- 
bination of measured quantities motivated by Eq. ([23)) . 
The common boost Ay t o of the quadrupole spectra is not 
observed for single-particle soft components — the mono- 
lithic boost distribution of the blast-wave model is in- 
consistent with data. The sharp edge of the quadrupole 
spectrum (narrow boost distribution) is particularly in- 



consistent with Hubble expansion, as shown in Fig. [12] 
Independence of the quadrupole component boost from 
the soft component is thus in conflict with model 1. 

E. Implications of the present analysis 

• Hadron y t spectra associated with the quadrupole 
component have been recovered from v^ipt) data 

• V2{pt) data trends are revealed as a complex inter- 
play of three hadron production mechanisms with 
accidental manifestations of mass dependence 

• The structure of the i>2{pt) ratio is dominated by 
the spectrum hard component above 0.5 GeV/c 

• Quadrupole hadrons come from a boosted source 
with narrow boost distribution not common to 
most hadrons; the number of hadrons from the 
quadrupole source is a small fraction of the total 

• The small-p t mass ordering invoked to support a 
hydro interpretation is a kinematic consequence of 
any common boosted source 

• The quadrupole component appears to be isolated 
from the rest of the collision evolution 

• Quadrupole spectra are substantially "cooler" than 
the single-particle spectrum soft component (i.e., 
the quadrupole and soft components are not in 
thermal equilibrium, with each other or with the 
spectrum hard component) 

• V2(j>t) trends interpreted as "constituent quark 
scaling" at intermediate pt do not relate to a hydro 
phenomenon or to hadron formation from a ther- 
malized partonic medium 

XVII. SUMMARY 

Elliptic flow (V2) measurements provide the primary 
support for claims of "perfect liquid" at RHIC. That cen- 
tral role motivates a careful re-examination of the inter- 
pretation of i>2 data in terms of hydrodynamic models. 
To that end I have reviewed azimuth correlation analy- 
sis methods and provided important generalizations. I 
described a method to extract quadrupole spectra on y t 
from i>2 (pt) data, and I used a limited V2 data sample for 
identified hadrons from minimum-bias Au-Au collisions 
at 200 GeV to illustrate properties of i>2 inferred from p t 
dependence and mass dependence. 

I reviewed an accurate two-component parameteriza- 
tion of hadron single-particle spectra on y t required to 
extract quadrupole spectra from V2(pt), and introduced 
the kinematics of boosted sources as an aid to interpret- 
ing features of «2 (pt ) ■ I expressed the functional form of 
V2{pt) as the product of two factors: p' t (pt in a boosted 
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frame) and the ratio of the sought-after quadrupole spec- 
trum to the single-particle spectrum 

I described the analysis steps required to combine the 
above elements so as to recover quadrupole spectra from 
Wa(pt) data and modeled the extracted quadrupole spec- 
tra with Levy distributions — boosted soft components 
S' NN (yt — Ay to). I compared the quadrupole spectrum 
component quantitatively to other spectrum components 
and to two hydro theory examples, and I estimated abso- 
lute quadrupole yields. Finally, I considered the impact 
of nonflow (minijet) contributions to v 2 measurements 

The conclusions from this analysis are as follows: 
Claims for v 2 scaling behavior supporting inference of a 
major role for constituent quarks in collision dynamics 
appear to be unsupported given the structure of ratio 
Vzipt) an d mixing of different physical mechanisms by 
that measure, especially above p t ~ 0.5 GeV/c. The true 
universality, as in Fig. [1] (right panel), is that of hadrons 
emitted from a common boosted source by the same 
hadronization mechanism as the single-particle spectrum 
soft component, albeit with a smaller "temperature." 
There is no support for a novel hadron production mech- 
anism. Monopole boost Ay t o is accurately obtained from 
Vzipt) data, but the (small) quadrupole absolute yield is 
inferred indirectly, since only the product of quadrupole 



boost Ay t 2 and absolute yield n c hi is measured directly. 

Analysis of data for three hadron species indicates that 
quadrupole yields relative to spectrum soft components 
are similar. The production mechanism for the soft- 
component yields in N-N collisions is the mechanism for 
the quadrupole yields in Au-Au collisions. Thus, only 
three numbers (two boosts and a ratio) are obtained from 
minimum-bias V^ijpt) data. Quadrupole (i^) data pro- 
vide no evidence for a thermalized system or for medium 
properties such as viscosity. The quadrupole component 
appears to result from an isolated dynamical process in- 
volving at most 5% of the hadrons in Au-Au collisions. 

The combination of those properties suggests that the 
azimuth quadrupole may be a new QCD phenomenon 
emerging at smaller QCD energy scales, interaction of 
QCD fields over large space-time volumes, which does not 
couple significantly to other collision processes and pro- 
duces a hadron spectrum significantly "cooler" than the 
spectrum soft component from nucleon fragmentation. 
The smaller slope parameter may result from reduced 
k t broadening of the QCD quadrupole field component 
compared to nucleon fragmentation (soft component). 

This work was supported in part by the Office of Sci- 
ence of the U.S. DoE under grant DE-FG03-97ER41020. 
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